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ABSTRACT 

We  study  the  properties  of  massive,  galactic- scale  outflows  of  molecular  gas  and  investigate  their  impact  on  galaxy  evolution.  We 
present  new  IRAM  PdBI  CO(l-O)  observations  of  local  ULIRGs  and  QSO  hosts:  clear  signature  of  massive  and  energetic  molecular 
outflows,  extending  on  kpc  scales,  is  found  in  the  CO(l-O)  kinematics  of  four  out  of  seven  sources,  with  measured  outflow  rates  of 
several  100  M0  yr_1.  We  combine  these  new  observations  with  data  from  the  literature,  and  explore  the  nature  and  origin  of  massive 
molecular  outflows  within  an  extended  sample  of  19  local  galaxies.  We  find  that  starburst-dominated  galaxies  have  an  outflow  rate 
comparable  to  their  SFR,  or  even  higher  by  a  factor  of  ~2-4,  implying  that  starbursts  can  indeed  be  effective  in  removing  cold  gas 
from  galaxies.  Nevertheless,  our  results  suggest  that  the  presence  of  an  AGN  can  boost  the  outflow  rate  by  a  large  factor,  which 
is  found  to  increase  with  the  LAGN/Lboi  ratio.  The  gas  depletion  time-scales  due  to  molecular  outflows  are  anti-correlated  with  the 
presence  and  luminosity  of  an  AGN  in  these  galaxies,  and  range  from  a  few  hundred  million  years  in  starburst  galaxies,  down  to 
just  a  few  million  years  in  galaxies  hosting  powerful  AGNs.  In  quasar  hosts  the  depletion  time-scales  due  to  the  outflow  are  much 
shorter  than  the  depletion  time-scales  due  to  star  formation.  We  estimate  the  outflow  kinetic  power  and  find  that,  for  galaxies  hosting 
powerful  AGNs,  it  corresponds  to  about  5%  of  the  AGN  luminosity,  as  expected  by  models  of  AGN  feedback.  Moreover,  we  find  that 
momentum  rates  of  about  20  LAGN/c  are  common  among  the  AGN-dominated  sources  in  our  sample.  For  “pure”  starburst  galaxies 
our  data  tentatively  support  models  in  which  outflows  are  mostly  momentum-driven  by  the  radiation  pressure  from  young  stars  onto 
dusty  clouds.  Overall,  our  results  indicate  that,  although  starbursts  are  effective  in  powering  massive  molecular  outflows,  the  presence 
of  an  AGN  may  strongly  enhance  such  outflows  and,  therefore,  have  a  profound  feedback  effect  on  the  evolution  of  galaxies,  by 
efficiently  removing  fuel  for  star  formation,  hence  quenching  star  formation. 

Key  words.  Galaxies:  active  -  Galaxies:  evolution  -  quasars:  general  -  Radio  lines:ISM  -  ISM:  molecules 


1.  Introduction 

The  recent  discovery  of  massive,  highly-energetic  and  kpc-scale 
molecular  outflows  in  local  AGN  hosts  and  ultra-luminous  in¬ 
frared  galaxies  (ULIRGs)  (Feruglio  et  al.  2010;  Fischer  et  al. 
2010;  Sturm  et  al.  2011;  Alatalo  et  al.  2011;  Aalto  et  al.  2012a; 


Cicone  et  al.  2012;  Dasyra  &  Combes  2012;  Nesvadba  et  al. 
2010,  2011;  Veilleux  et  al.  2013;  Spoon  et  al.  2013;  Gonzalez- 
Alfonso  et  al.  2013)  has  provided  important  evidence  for  neg¬ 
ative  feedback  on  star  formation  in  action  in  galaxies  (see  also 
the  review  by  Fabian  2012).  The  major  breakthrough  achieved 
by  these  studies  is  the  confirmation  that,  whatever  is  the  driving 
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mechanism  (AGN  or  star  formation  itself),  these  outflows  affect 
the  phase  of  the  interstellar  medium  (ISM)  out  of  which  stars 
forms,  i.e.  the  cold  molecular  gas.  Therefore,  such  feedback 
mechanism  can  have  a  significant  impact  on  the  evolution  of  the 
host  galaxy,  by  regulating,  through  outflows,  the  amount  of  cold 
and  dense  gas  available  in  the  ISM  for  star  formation  and  black 
hole  accretion.  In  particular,  negative  feedback  from  quasars  is 
believed  to  prevent  massive  galaxies  from  overgrowing  (hence 
explaining  the  steep  decline  of  the  stellar  mass  function  at  high 
masses),  to  account  for  the  “red- and- dead”  properties  of  massive 
ellipticals  (Cowie  et  al.|1996[|Baldry  et  al.|2004[|Perez-Gonzalez| 


|et  al.|2008),  and  to  drive  the  correlations  between  supermassive 

black  holes  and  bulge  properties  that  are  observed  in  the  local 
Universe  (e.g.  Magorrian  et  al.||1998|  |Marconi  &  Hunt||2003 


Ferrarese  &  Ford|2005[|Kormendy  &  Ho|2013|). 


Indeed,  most  theoretical  models  propose  that  the  super  mas¬ 
sive  black  hole,  during  its  bright  active  phase  (“quasar-mode" 
feedback),  characterised  by  a  high  accretion  rate,  can  expel  the 
ISM  out  of  the  host  galaxy,  and  eventually  clears  the  galaxy  of  its 
cold  gas  reservoir,  although  the  details  of  the  mechanism  respon¬ 
sible  for  coupling  the  quasar  radiation  with  the  galactic  ISM  are 


still  debated  (e.g.  Silk  &  Rees  1998;  Granato  et  al.  2004;  Lapi 

et  al. 

2005 

Di  Matteo  et  al.  2005 ;  Springel  et  al.  2005 ;  Croton 

et  al. 

2006 

Hopkins  et  al. 

2008; 

Menci  et  al. 

2008;  Narayanan 

et  al. 

2008;  Fabian  et  al. 

2009; 

King  2010f 

Zubovas  &  King 

2012 

Faucher-Giguere  &  Quataert  20 1 2 ). 

While  various  techniques  have  been  used  to  identify  other 
phases  of  outflows  in  active  galaxies,  both  locally  (e.g.  |Rupke 


&  Veilleux  2011, 20 13b  [ 

Muller-Sanchez  et  al.  2011;  Mullaney 

et  al.  2013,  Emonts  et  al. 

2005;  Morganti  et  al.  2003,  2005  ;  Ro- 

driguez  Zaurin  et  al.  2013)  and  at  high  redshift  (e.g.  Harrison 

et  al.  2012;tCano-Diaz  et 

al.  2012;  Diamond- Stanic  et  al.  2012; 

Bradshaw  et  al.  112013]  [Alexander  et  al.||2010[  |Maiolino  et  al.| 

2012J,  the  interferometric  mapping  of  the  high  velocity  compo¬ 
nent  of  CO  millimeter  emission  has  proved  to  be  very  effective. 
Indeed  this  method,  as  mentioned,  allows  us  to  directly  trace  out¬ 
flows  of  the  molecular  phase  out  of  which  stars  form,  which  in 
most  galaxies  also  represents  the  bulk  of  the  ISM;  moreover  the 
spectro-imaging  information  delivered  by  millimeter  interferom¬ 
eters  allows  us  to  determine  the  size  and  geometry  of  the  outflow, 
hence  to  accurately  estimate  the  outflow  rate  and  energetics. 

In  this  work  we  aim  to  address  some  open  questions,  such 
as:  are  massive  molecular  outflows  effective  in  quenching  star 
formation  in  the  host  galaxy?  Is  the  central  AGN  the  domi¬ 
nant  power  source  and,  consequently,  are  AGNs  controlling  the 
growth  and  evolution  of  galaxies?  Are  current  AGN-feedback 
models  adequate  to  explain  observations? 

In  this  paper  we  present  new  IRAM  PdBI  CO(l-O)  observa¬ 
tions  for  a  sample  of  7  local  active  galaxies,  aimed  at  identifying 
the  presence  of  molecular  outflows.  Our  new  data  are  combined 
with  previous  results  from  our  and  other  teams,  with  the  goal  of 
characterising  the  physical  properties,  origin  and  nature  of  mas¬ 
sive  molecular  outflows,  and  their  impact  on  galaxy  evolution. 

A  //0=70.4  km  s-1  Mpc-1,  =  0.27,  Da  =  0.73  cosmol¬ 
ogy  is  adopted  throughout  this  work. 


2.  Sample  Selection  and  IRAM  PdBI  Observations 

We  observed  seven  local  active  luminous  and  ultra-luminous  in¬ 
frared  galaxies  in  their  CO(l-O)  molecular  transition  with  the 
IRAM  Plateau  de  Bure  Interferometer.  Since  these  are  still  ex¬ 
ploratory  observations,  the  sample  has  a  mixed  composition. 
Out  of  the  seven  sources,  five  are  extracted  from  the  SHIN¬ 
ING  project,  which  is  a  Herschel-PACS  spectroscopic  guar¬ 


anteed  time  key  program,  targeting  starbursts,  Seyfert  galax¬ 
ies,  and  (ultra)  luminous  infrared  galaxies  ([Fischer  et  al.|2010"| 
|Sturm  etal.|2011H  Veilleux  et  al.|2013|).  Four  of  these  ULIRGs, 
namely  Mrk  273,  IRAS  F08572+3915,  IRAS  F10565+2448  and 
IRAS  23365+3604,  show  evidence,  in  their  far-IR  spectrum,  for 
P-Cygni  profiles  of  the  OH  transitions  at  79  and  119yum,  which 
unambiguously  trace  molecular  outflows  (Sturm  et  al.||201l] 
|Veilleux  et  al.|2013|).  In  the  other  SHINING  source,  Mrk  876,  as 
well  as  in  the  quasar  I  Zw  1 ,  also  observed  with  PACS  in  a  cycle 
1  open-time  program,  OH  is  detected  only  in  emission  (Veilleux 
|et  al.|2013|.  These  two  sources  are  incorporated  in  our  sample  to 
investigate  whether  evidence  of  outflow  is  observed  in  CO  emis¬ 
sion,  hence  to  test  whether  the  detection  of  OH  blue-shifted  ab¬ 
sorption  is  a  pre-requisite  for  the  presence  of  molecular  outflows. 
Finally,  IRAS  F23060+0505  is  one  of  the  most  powerful  type  2 
quasars,  and,  by  observing  it,  we  aim  to  investigate  the  presence 
of  molecular  outflows  in  this  class  of  very  powerful,  but  heavily 
embedded  objects.  The  list  of  targets  is  given  in  Table [T] 

Observations  of  the  CO(l-O)  rotational  transition  (vrest  = 
115.271  GHz)  were  performed  with  the  IRAM  Plateau  de  Bure 
Interferometer  (PdBI)  between  May,  2010  and  July,  2012.  A 
technical  description  of  the  observational  parameters  is  provided 
in  Table  [T]  In  the  compact  (D)  array  configuration,  only  5  of 
the  15m  antennas  were  used,  while  observations  in  C  config¬ 
uration  were  carried  out  with  the  complete  array  (6  antennas). 
Data  calibration  and  analysis  were  performed  by  using  the  CLIC 
and  MAPPING  softwares  within  the  GILDAS  package.  The  flux 
calibration  accuracy  at  these  frequencies  is  about  10%  ([Castro- 


|Carrizo  &  Neri|20l0[).  We  exploit  the  wide-band  (WideX)  cor¬ 
relator  offered  by  the  PdBI,  which  provides  a  spectral  resolution 
of  1.95  MHz  (corresponding,  at  an  average  observed  frequency 
of  — 1 10  GHz,  to  channel  widths  of  ~  5.3  km  s_1)  over  its  full 
bandwidth  of  3.6  GHz.  For  display  purposes,  we  bin  the  spectra 
using  a  re-binning  factor  that  varies,  from  source  to  source,  from 
5  to  20  channels. 


3.  Results 

3.1.  Criterion  for  Molecular  Outflow  Detection 

The  adopted  criterion  to  claim  the  detection  of  a  molecular  out¬ 
flow  is  that  at  least  two  of  the  following  conditions  must  be  sat¬ 
isfied: 

1.  CO(l-O)  emission  with  velocities  higher  than  500  km  s-1  is 
detected  in  the  interferometric  continuum- subtracted  maps 
with  a  significance  of  at  least  5<x; 

2.  Broad  CO(l-O)  wings  with  velocities  higher  than  300  km 
s-1  are  detected  in  the  interferometric  continuum- subtracted 
maps  with  a  significance  of  at  least  5<x  and  are  found  to  de¬ 
viate  from  the  rotational  pattern  (for  example  in  the  position- 
velocity  diagram); 

3.  A  molecular  outflow  has  been  already  detected  through  P- 
Cygni  profiles  of  molecular  transitions  (e.g.  OH,  H2O). 

This  criterion  selects  only  those  galaxies  for  which  unambiguous 
evidence  for  molecular  outflows  is  found;  sources  which  only 
meet  conditions  1  or  2  may  also  be  hosting  a  powerful  molecular 
outflow,  but  additional  observational  evidence  is  needed  in  order 
to  exclude  that  the  high  velocity  CO  emission  is  tracing  other 
mechanisms,  such  as  gas  rotation,  infalls  or  also  disturbed  gas 
kinematics  due  to  an  ongoing  galaxy  merging. 

According  to  the  criterion  outlined  above,  we  find  evidence 
for  massive  molecular  outflows  in  the  CO(l-O)  data  for  four  out 
of  seven  galaxies  of  our  sample,  namely:  IRAS  F08572+3915, 
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Table  1:  Description  of  the  IRAM  PdBI  CO(l-O)  Observations 


Object 

RAO 

(J2000) 

DECO 

(J2000) 

Date(s) 

Obs.  Freq. 
(GHz) 

Min,  Max 
Baseline  (m) 

Conf. 

On  Source 
Time  (hrs) 

FoV 

(arcsec) 

Synt.  Beam 
(arcsec) 

IRAS  F08572+3915 

09:00:25.37 

39:03:53.68 

May-Oct  2011 

108.930 

(15.0,  176.0) 

C+D 

20 

46.3  x  46.3 

3.1  x  2.7 

IRAS  F10565+2448 

10:59:18.14 

24:32:34.42 

Jun-Oct  2011 

110.509 

(16.5,  175.2) 

C+D 

10 

45.6  x  45.6 

3.5  x  2.9 

IRAS  23365+3604 

23:39:01.23 

36:21:09.30 

Jun-Oct  2011 

108.299 

(17.0,  175.8) 

C+D 

10 

46.5  x  46.5 

3.1  x  2.6 

Mrk  273 

13:44:42.16 

55:53:13.21 

Jun-Jul  2012 

111.076 

(17.5,  102.5) 

D 

12 

45.4  x  45.4 

5.0  x  4.0 

IRAS  F23060+0505 

23:08:33.93 

05:21:29.90 

May  2010 

98.270 

(15.0,  94.5) 

D 

3 

51.3x51.3 

5.3  x  4.5 

Mrk  876 

16:13:57.21 

65:43:10.60 

May  2010 

102.100 

(16.6,  93.2) 

D 

6 

49.4  x  49.4 

6.0  x  3.9 

I  Zw  1 

00:53:34.94 

12:41:36.20 

May  2010 

108.629 

(17.8,  95.6) 

D 

8 

46.4  x  46.4 

4.7  x  3.7 

IRAS  F10565+2448,  IRAS  23365+3604  and  Mrk  273,  which 
also  exhibit  prominent  OH  (79  and  119yum)  P-Cygni  profiles  in 
their  Herschel-PACS  spectra  ([Sturm  et  al.|2011[  |Veilleux  et  al.| 
|2013|).  This  proves  that,  once  that  the  presence  of  a  molecular 
outflow  has  been  assessed  by  the  detection  of  blue-shifted  OH 
absorption,  interferometric  CO  observations  are  a  very  success¬ 
ful  complementary  technique  to  study  more  in  detail  the  char¬ 
acteristics  of  the  molecular  outflow,  allowing  us  to  directly  esti¬ 
mate  the  mass  and  size  of  the  outflow  and  hence  derive  other  im¬ 
portant  properties  (e.g.  outflow  rates,  kinetic  energy,  momentum 
rates  etc.).  In  the  powerful  type  2  quasar  IRAS  F23060+0505 
we  detect,  at  high  significance,  broad  CO(l-O)  wings  with 
|v|  >  300  km  s-1,  which  deviate  from  the  rotation  pattern  in  the 
PV  diagram.  IRAS  F23060+0505  therefore  meets  one  of  the 
conditions  mentioned  above  and  is  very  likely  to  host  a  mas¬ 
sive  molecular  outflow.  However,  since  gas  at  |v|  >  500  km  s-1 
is  only  marginally  detected,  and  since  IRAS  F23060+0505  was 
not  observed  by  Herschel,  this  source  does  not  fully  qualify  as 
a  molecular  outflow  “detection”.  The  remaining  two  AGN-host 
galaxies  of  our  sample,  Mrk  876  and  I  Zw  1,  do  not  fully  meet 
any  of  the  previously  mentioned  conditions  for  the  detection  of 
molecular  outflows,  although  in  the  case  of  Mrk  876  we  do  ob¬ 
serve  some  marginal  evidence  for  high  velocity  and  spatially  ex¬ 
tended  CO(l-O)  emission.  In  the  absence  of  additional  observa¬ 
tional  evidence  and  of  better  quality  data,  we  use  our  CO(l-O) 
observations  to  estimate,  for  IRAS  F23060+0505,  Mrk  876  and 
I  Zw  1,  upper  limits  on  the  mass  of  molecular  outflowing  gas,  on 
the  outflow  rates  and  other  outflow  properties. 


3.2.  Deriving  Physical  Properties  of  Molecular  Outflows 


Our  approach  to  the  analysis  of  the  broad  wings  of  the  CO(l- 
0)  emission  lines  in  the  seven  sources  is  similar  to  that  adopted 
for  studying  the  massive  outflows  detected  in  the  CO(l-O)  and 
CO(2-l)  transitions  in  Mrk  231  ([Feruglio  et  al.||2010t  |Cicone| 
et  al.||2Q12|)  and  in  the  [CII]  158//m  line  in  the  high  redshift 


quasar  SDSS  J1 148+5 152  ([Maiolino  et  al.|2012|). 


Our  method  for  estimating  the  integrated  flux  and  the  spa¬ 
tial  extension  of  the  CO(l-O)  wings  relies  almost  exclusively  on 
the  analysis  of  the  uv  visibility  data,  therefore  producing  results 
which  are  independent  of  aperture  effects  and  of  the  cleaning 
process.  More  specifically,  we  first  estimate  the  3mm  continuum 
emission  from  the  line-free  frequency  ranges  in  our  wide-band 
(WideX)  observations,  and  then  subtract  it  directly  from  the  uv 
data.  We  then  average  the  frequency  channels  corresponding  to 
the  CO(l-O)  wings,  therefore  producing  uv  tables  for  both  the 
blue  and  the  red- shifted  CO  wing,  which  are  used  to  estimate 
the  integrated  flux  and  the  size  of  the  CO  wings.  To  infer  the 
flux  (both  for  the  core  of  the  line  and  for  the  wings),  we  bin  the 
visibility  vs  uv  radius  plots  (form  here  simply  referred  to  as  “wv” 
plots)  in  baseline  steps  of  20m,  and  use  the  visibility  amplitude 
of  the  first  point  of  this  plot,  which  approximates  the  total  flux 
(i.e.  the  flux  that  would  be  measured  at  a  zero  antenna  separa¬ 


tion).  To  infer  the  spatial  extent  of  the  wings,  we  fit  the  uv  plot 
with  a  circular  Gaussian  model,  and  use  the  resulting  FWHM  as 
an  estimate  of  the  size  (diameter)  of  the  CO  wings  emission.  We 
stress  that  this  technique,  based  on  the  uv  data  rather  than  on  the 
interferometric  deconvolved  maps,  allows  us  to  estimate  the  size 
of  structures  which  are  comparable  to  the  synthesised  beam,  on 
condition  that  the  signal  from  the  source  is  detected  at  least  at  a 
8 cr  significance.  Moreover,  the  advantage  of  working  directly  in 
the  uv  plane  for  estimating  the  fluxes  and  the  sizes,  is  that  at  this 
level  no  alterations  are  made  on  the  flux  distribution  as  conse¬ 
quence  of  the  “cleaning”  process.  The  fluxes  and  sizes  inferred 
for  the  broad  CO  wings  of  our  seven  sources  are  listed  in  Table 
[5]  For  I  Zw  1,  where  CO(l-O)  broad  wings  are  not  detected^, 
we  provide  upper  limits. 

We  use  the  quantities  listed  in  Table  [2]  to  derive  physi¬ 
cally  relevant  outflow  properties:  outflow  molecular  gas  masses, 
mass-loss  rates,  molecular  gas  depletion  time-scales  due  to  the 
outflow,  kinetic  powers  and  momentum  rates  of  the  molecular 
outflows,  as  detailed  in  the  following.  Outflow  rates  and  ener¬ 
getics  are  estimated  by  adopting  the  same  method  as  for  Mrk 
231  ([Feruglio  et  aH2010|,  NGC  6240  ([Feruglio  et  al.[2013a|b|) 


and  SDSS  J1 148+5152  ([Maiolino  et  al.l2012[).  We  refer  to  those 
papers  for  an  exhaustive  description  of  our  outflow  model  and  its 
implications.  Here,  we  only  briefly  summarise  the  method. 

The  excitation  study  of  the  massive  molecular  outflow  in  the 
ULIRG  Mrk  23 1  indicated  that  the  bulk  of  the  molecular  gas  in 
the  galaxy  and  the  outflowing  gas  have  roughly  similar  excitation 
([Cicone  etal.|2012|  Feruglio  et  al.  in  prep.).  A  direct  implication 
of  this  result  is  that,  to  evaluate  the  molecular  outflow  mass  from 
the  CO  luminosity  of  the  broad  wings,  a  better  choice  would  be 
to  use  the  standard  ULIRG  CO-to-H2  conversion  factor,  rather 


than  the  lower,  more  conservative  value  assumed  in  Feruglio 
[et  al.|([20T()[).  A  similar  result  has  been  found  in  M82,  where  a 
detailed  modelling  of  the  conversion  factor  of  the  gas  in  the  out- 
11  ow  has  been  performed  by|Weifi  et  al.|(|200T])  and | Walter  et  al. 
(2002).  As  a  consequence,  for  all  molecular  outflows  treated  in 
this  paper,  we  adopt  trco(i-o)=  0-8  M0(Kkms_1  pc1 2)-1. 


As  already  mentioned,  for  the  galaxies  of  our  PdBI  sam¬ 
ple,  as  well  as  for  Mrk  231  and  NGC  6240,  we  assume  the 
same  spherical,  or  multi-conical,  geometry  as  in|Maiolino  et  al. 
(|2012[).  As  discussed  in  the  latter  paper,  if  the  outflowing  clouds 
are  assumed  to  uniformly  populate  the  spherical  (or  multi- 
conical)  region  affected  by  the  outflow,  then  the  outflow  rate  is 
given  by  the  relation: 


Mqf  =  3v  (Mqf/Rof)  (1) 


1  We  note  that  throughout  this  paper,  the  equivalence  “broad  CO  wings 

=  molecular  outflow”  is  not  valid,  because  of  the  criterion  that  we  define 
for  the  molecular  outflow  detection  (Section  3.1).  More  specifically,  al¬ 
though  IRAS  F23060+0505,  Mrk  876,  and  I  Zw  1  are  all  treated  as  non 
detections  for  the  molecular  outflow,  IRAS  F23060+0505  and  Mrk  876 
differ  from  I  Zw  1  in  that  broad  CO  wings  are  (marginally)  detected, 
and  we  can  estimate  their  flux  and  size  from  our  CO(l-O)  observations. 
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The  formula  above  provides  an  appropriate  description  if  the  out¬ 
flow  is  continuously  refilled  with  clouds  ejected  from  the  galac¬ 
tic  gaseous  disc,  as  expected  in  feedback  scenarios.  However,  if 
the  outflow  is  instead  associated  with  a  single  explosive  event, 
in  which  the  clouds  are  ejected,  then  the  outflow  rate  is  more 
properly  described  by  the  gas  mass  of  the  clouds  divided  by 
the  dynamical  time  required  by  the  clouds  to  reach  the  current 
location,  resulting  in  an  outflow  rate  which  is  one  third  of  the 
one  given  by  Eq.  [T]  Higher  angular  resolution  data  would  be 
required  to  discriminate  between  these  two  possible  scenarios, 
and  in  particular  to  clearly  distinguish  between  a  uniform  dis¬ 
tribution  of  clouds  within  outflow  cones  (former  scenario,  de¬ 
scribed  by  Eq.[T]),  or  a  shell-like  geometry  (latter  scenario).  The 
data  available  so  far  for  our  sources  (including  the  well- studied 
Mrk  231),  are  more  consistent  with  a  uniform  outflow,  rather 
than  with  the  shell-like  explosive  scenario,  hence  in  this  work 
we  adopt  the  prescription  given  by  Eq.  [T]  However,  the  reader 
should  be  aware  that  if  future  data  favour  an  explosive-like  (i.e. 
shell-like)  scenario,  then  the  figures  for  outflow  rates,  kinetic 
powers  and  momentum  rates  may  need  to  be  revised  downward. 
For  some  objects  taken  from  the  literature  and  discussed  in  Sec¬ 
tion  4  this  seems  indeed  to  be  the  case  (and  this  is  properly  taken 
into  account  for  the  determination  of  the  outflow  rate  for  these 
sources).  In  Eq.  [T]  we  conservatively  use  the  average  velocity 
in  the  range  used  to  integrate  the  broad  CO  wings  (rather  than 
the  maximum  velocity,  both  listed  in  Table  [4]).  We  note  that,  if 
the  outflowing  gas  intersects  our  line  of  sight  (which  is  always 
true,  as  probed  by  the  Herschel-PACS  detection  of  OH  P-Cygni 
profiles),  Eq.  [T]  would  provide  a  fairly  accurate  estimate  of  the 
outflow  rate  also  in  the  case  of  a  bi-conical  wind  geometry.  The 
kinetic  power  and  momentum  rate  of  the  outflow  are  simply  cal¬ 
culated  as  1  /2  v2 *  Mh2,  of  and  v  Mh2,  of,  respectively. 

All  the  outflow  properties  described  above  are  listed  in  Table 
[4]  for  the  extended  sample  of  galaxies,  consisting  of  our  seven 
sources  observed  with  the  PdBI,  and  of  an  heterogeneous  sample 
of  12  local  galaxies  hosting  molecular  outflows,  drawn  from  the 
literature  (further  explanation  in  Section  4). 

For  IRAS  F23060+0505  and  Mrk  876,  we  use  the  values  re¬ 
ported  in  Table [2]for  flux  and  size  of  the  CO  wings  to  infer  upper 
limits  on  the  outflow  properties.  In  the  case  of  I  Zw  1,  the  up¬ 
per  limits  on  the  outflow  parameters  are  calculated  by  using  the 
upper  limit  that  we  estimate  for  the  flux  in  the  CO  wings,  listed 
in  Table [2|  and  hypothetical  values  for  the  outflow  extension  (ra¬ 
dius)  and  velocity  of  0.5  kpc  and  500  km  s-1,  respectively. 

3.3.  Analysis  of  the  I  RAM  PdBI  CO(I-O)  Observations 

In  the  following  we  present  our  IRAM  PdBI  CO(l-O)  observa¬ 
tions  of  the  seven  local  AGN-host  galaxies  listed  in  Table  [T 
mainly  focusing  on  the  evidence  for  broad  CO(l-O)  wings  anc 
their  properties.  For  additional  information  about  the  individual 
sources  and  further  discussion  of  our  CO(l-O)  observations  we 
refer  to  Appendix  A. 


IRAS  F08572+3915 

The  continuum- subtracted  CO(l-O)  emission  line  profile  of 
IRAS  F08572+3915  is  shown  in  Figure  [lj  in  which  we  report 
the  PdBI  spectra  extracted  from  an  aperture  of  diameter  ^  10 
arcsec  (a)  and  from  the  central  beam  (b),  which  is  actually  sim¬ 
ply  measured  from  the  central  pixel  (since  the  maps  are  produced 
in  units  of  mJy/beam).  In  this  source  we  estimated  the  contin¬ 
uum  emission  over  the  velocity  ranges  ve(-3700,  -2500)  km  s-1 


and  ve  (2500,  4000)  km  s-1.  We  measure  a  total  integrated  line 
flux  of  S co, tot  =  (10.80  ±  0.70)  Jy  km  s-1  within  ±  2000  km 
s-1  from  the  systemic  velocity,  which  is  consistent  with  previous 
single  dish  and  interferometric  observations  ([Evans  et  al.|2002 ; 


Solomon  et  al.[1997t|Papadopoulos  et  al.|2012b|). 


The  large  bandwidth  of  WideX  allows  us  to  detect,  for  the 
first  time  in  this  source,  spectacular  broad  wings  of  the  CO(l-O) 
emission  line,  tracing  an  outflow  with  velocities  of  up  to  1200 
km  s-1.  The  two  spectral  profiles  presented  in  Figure  0a-b) 
show  interesting  differences.  In  the  spectrum  extracted  from  the 
10  arcsec  aperture,  the  blue  and  the  red  wing  appear  almost  per¬ 
fectly  symmetric  with  respect  to  the  line  core,  and  their  emis¬ 
sion  peaks  at  velocities  of  about  ±900  km  s-1  from  the  systemic. 
Moreover,  the  narrow  core  can  be  fitted  quite  well  with  a  sin¬ 
gle  Gaussian  function.  Conversely,  in  the  spectrum  extracted 
from  the  central  pixel,  the  fit  to  the  narrow  core  shows  a  clear 
excess  of  emission  at  ve  (-200,  -400)  km  s-1,  which  is  not  dis¬ 
cernible  in  the  10  arcsec  aperture  spectrum,  probably  because 
of  the  lower  signal-to-noise.  Nonetheless,  the  most  obvious  dif¬ 
ference  between  the  two  spectra  is  the  appearance  of  the  high 
velocity  wings:  in  the  central  pixel  spectrum,  while  the  blue- 
shifted  wing  appears  very  prominent,  the  red-shifted  emission  at 
v  >  500  kms-1  is  not  recovered,  which  would  suggest  that  this 
component  is  mostly  extended.  We  note  that  some  very  recent 
follow-up  observations  of  this  source,  obtained  by  us  in  PdBI 
A-configuration,  resolve  out  the  red  wing,  hence  confirming  that 
this  must  trace  a  very  extended  outflow. 


The  cleaned  maps  of  the  blue  and  red-shifted  wings  are  pre¬ 
sented  in  panels  (c-d)  of  Fig.  [T]  Both  the  wings  are  detected  at  a 
>  10cr  significance  in  the  maps,  which  is  also  the  significance 
of  detections  of  wings  at  |v|  >  5OOkms_10.  The  molecular 
outflow  is  detected  in  correspondence  of  the  north-west  galaxy 
of  the  merging  system,  confirming  findings  of  previous  studies. 
The  peak  of  the  red-shifted  emission  is  slightly  offset  (~  1.4  arc¬ 
sec)  with  respect  to  the  pointing  and  phase  centre  (indicated  by 
the  cross).  Moreover,  the  red  wing  map  exhibits  a  structure  ex¬ 
tended  to  the  north,  which  further  confirms  the  hypothesis  of  a 
very  extended  red-shifted  emission.  The  fluxes  measured  in  the 
blue  and  red-shifted  wings  of  the  CO(l-O)  emission  line  are  re¬ 
ported  in  Table  [2]  The  analysis  of  the  uv  data  gives  an  average 
size  (FWHM)  of  1.56  ±  0.33  kpc  for  the  massive  molecular  out¬ 
flow  in  IRAS  F08572+3915  (see  Appendix  A  for  the  details). 
The  velocity  map  of  the  central  core  emission  (shown  in  the  Ap¬ 
pendix  A)  shows  a  regular  rotation  pattern,  with  a  major  axis  ori¬ 
ented  east-west,  similarly  to  the  narrow  ionised  gas  component 
(|Rupke  &  Veilleux|2013b|).  The  position- velocity  diagram  along 
the  rotation  major  axis  is  presented  in  panels  (e-f)  of  Fig.  [T| 
The  inset  (f)  shows  the  central  rotation  pattern,  which  is  instead 
burned  out  in  panel  (e)  with  the  adopted  colour  cuts.  The  signa¬ 
ture  of  the  extreme  molecular  outflow  is  clearly  seen  through  the 
high  velocity  gas  (along  the  horizontal  axis  of  the  plot)  that  does 
not  follow  the  rotation  curve  of  the  galaxy. 


Summarising,  the  analysis  of  our  PdBI  CO(l-O)  observations 
of  the  Sy2-UFIRG  IRAS  F08572+3915  convincingly  reveals  the 
presence  of  a  powerful  and  massive  molecular  outflow,  which 
extends  on  both  sides  of  the  galaxy  up  to  velocities  of  -1200 
km  s-1  and  on  physical  scales  of  -1  kpc,  carrying  a  mass  of 
molecular  gas  of  about  4  x  108  M0,  resulting  in  an  outflow  mass- 
loss  rate  of  -  1200  M©  yr-1. 


2  By  restricting  the  velocity  range  of  integration  for  the  blue  wing  to 

velocities  v  <  -500  km  s-1  we  obtain,  similarly,  a  S/N=12  detection. 
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(a)  IRAS  FG8572+3915  (10  arcsec)  (b)  IRAS  F08572+3915  (cen  pix) 


Fig.  1:  Continuum-subtracted  IRAM-PdBI  spectra,  maps  and  position- velocity  diagram  of  the  CO(l-O)  emission  line  of  IRAS 
F08572+3915.  Panels  (a)  and  (b)  show  the  spectra  extracted  from  a  circular  aperture  with  diameter  of  10  arcsec  and  from  the 
central  pixel  (i.e.  the  centroid  of  the  integrated  CO  emission),  respectively.  For  display  purposes,  the  spectra  are  re-binned  in 
channels  of  107  km  s-1  (a)  and  54  km  s-1  (b).  In  both  panels  the  narrow  core  of  the  line  is  fitted  with  a  single  Gaussian  function 
(red  dashed  line).  Note  the  clear  detection  of  broad  wings  in  the  spectra,  (c-d)  Map  of  the  CO(l-O)  emission  integrated  in  the  blue 
and  red  wings.  Contours  correspond  to  ler  (l<x  rms  level  of  the  two  maps  is  0.1  mJy  beam-1).  The  physical  scale  at  the  redshift  of 
the  source  is  1.122  kpc  arcsec-1.  The  cross  indicates  the  peak  of  the  radio  (VLB I)  emission,  which  corresponds  to  the  NW  nucleus 
of  the  merging  pair  IRAS  F08572+3915.  (e-f)  Position-velocity  diagram  along  the  major  axis  of  the  molecular  disk  rotation  (as 
traced  by  the  narrow  core  of  the  CO(l-O)  line).  The  inset  (f)  highlights  the  rotation  pattern  traced  by  the  narrow  component  of  the 
CO(l-O)  line.  Contours  are  in  steps  of  l<x  (up  to  10<x)  in  panel  (e)  and  5 cr  in  panel  (f);  negative  contours  are  in  steps  of  1  cr  (magenta 
dashed  lines).  Each  pixel  corresponds  to  0.9  arcsec  x  90  km  s-1  in  diagram  (e),  and  to  0.3  arcsec  x  30  km  s-1  in  the  inset  (f).  Note 
the  presence  of  high  velocity  gas,  which  does  not  follow  the  rotation  pattern  and  extends  to  velocities  as  high  as  1200  km  s-1.  This 
traces  the  same  molecular  outflow  detected  in  OH  by  Herschel-PACS. 


IRAS  FI 0565+2448 

The  previous  CO(l-O)  observations  of  the  Sy2-ULIRG  IRAS 
FI 05 65 +244 8,  due  to  their  narrow  bandwidth,  could  not  detect 
the  broad  component  superposed  on  the  narrow  core  of  the  line. 
This  broad  component  is  clearly  visible  in  our  new  PdBI  spectra 
shown  in  Figure  (2ja-b).  The  broad  wings  in  this  source  are  de¬ 
tected  up  to  about  v  =  +  600  km  s-1.  The  continuum  emission 
was  estimated  in  the  velocity  ranges  ve(-3600,  -1200)  km  s-1 
and  ve  (1200,  4000)  km  s-1.  We  stress  that  the  choice  of  a  dif¬ 
ferent  continuum  window  with  respect  to  IRAS  F08572+3915 
and  to  the  other  sources,  does  not  introduce  artificial  bias  on  the 
shape  and  velocity  extension  of  the  broad  wings.  Indeed  for  each 
source  we  first  inspect  the  appearance  and  shape  of  the  broad 
wings  in  the  non  continuum- subtracted  spectrum. 

We  note  that  the  spectrum  extracted  from  a  circular  aper¬ 
ture  with  a  diameter  of  10  arcsec  (Fig  |2^)  and  the  one  ex¬ 
tracted  from  the  central  pixel  (Fig [2}))  show  similar  line  profiles. 
We  measure  a  total  integrated  flux  of  Sco,tot  =  (108.0  +  1.0) 


Jy  km  s-1,  which  is  significantly  (>  50  %)  larger  than  the  flux 
recovered  by  Downes  &  Solomon  (1998)  and  ~17%  higher  than 
the  value  given  by|Chung  et  al.|{|2009 )  (see  also  [Papadopoulos 
|et  al.|(|2012b|)).  We  report  in  Fig.  |2[c-d)  the  blue  and  red  wing 
maps  obtained  by  integrating  the  CO(l-O)  emission  in  the  ve¬ 
locity  ranges  (-600,  -300)  km  s-1  and  (300,  600)  km  s-1.  CO 
wings  at  |v|  >  300  km  s-1  are  detected  with  a  S/N  =12,  and  the 
position- velocity  diagram  along  the  major  axis  of  rotation  (pan¬ 
els  (e-f)  of  Fig.  [2]),  proves  that  this  high  velocity  gas  does  deviate 
from  the  central  rotation  pattern.  Because  Herschel-PACS  ob¬ 
servations  revealed,  in  this  source,  prominent  OH  P-Cygni  pro- 
files  with  blue-shifted  velocities  of  up  to  950  km  s-1  (Veilleux 


et  al.||2013]),  the  high  significance  detection  of  CO  wings  with 
|v|  >  300  km  s-1  deviating  from  the  rotational  pattern  is  suffi¬ 
cient,  according  to  our  criterion  (Section  3.1),  to  claim  the  de¬ 
tection  of  a  massive  molecular  outflow.  We  note  that  very  high 
velocity  (|v|  >  500  km  s-1)  CO(l-O)  emission  is  detected  in  this 
galaxy  at  a  significance  level  of  4 cr. 
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Fig.  2:  Continuum-subtracted  IRAM-PdBI  spectra,  maps  and  position- velocity  diagram  of  the  CO(l-O)  emission  line  of  IRAS 
F10565+2448.  Panels  (a)  and  (b)  show  the  spectra  extracted  from  a  circular  aperture  with  diameter  of  10  arcsec  and  from  the 
central  pixel,  respectively.  For  display  purposes,  the  spectra  are  re-binned  in  channels  of  27  km  s-1.  In  both  panels  the  narrow  core 
of  the  line  is  fitted  with  a  single  Gaussian  function  (red  dashed  line).  Note  the  appearance,  in  both  spectra,  of  broad  CO  wings  up  to 
±600  km  s-1.  (c-d)  Cleaned  map  of  the  emission  integrated  in  the  blue  and  red-shifted  CO(l-O)  wings.  Contours  correspond  to  lcr 
(1  cr  rms  level  of  the  two  maps  is  0.2  mJy  beam-1).  The  cross  indicates  the  peak  of  the  radio  (VLBI)  emission.  The  physical  scale 
at  the  redshift  of  the  source  is  0.846  kpc  arcsec-1.  (e-f)  Position-velocity  diagram  along  the  major  axis  of  rotation.  The  inset  (f) 
highlights  the  rotation  pattern  traced  by  the  narrow  component  of  the  CO(l-O)  line.  Contours  are  in  steps  of  lcr  (up  to  10cr)  in  panel 
(e),  and  5cr  in  panel  (f)  (starting  from  10cr  in  panel  (f));  negative  contours  are  in  steps  of  lcr  (magenta  dashed  lines).  Each  pixel 
corresponds  to  1.3  arcsec  x  53  km  s-1  in  diagram  (e),  and  to  0.6  arcsec  x  26  km  s-1  in  the  inset  (f).  Note  the  clear  signature  of  gas 
at  velocities  |v|  >  300  km  s-1,  which  does  not  follow  the  rotation  curve  and  is  therefore  ascribed  to  the  molecular  outflow  detected 
by  Herschel-PACS. 


Our  observations  resolve  the  outflow  traced  by  the  broad 
wings  of  the  CO(l-O)  line,  and  the  fit  to  the  uv  data  results  in 
an  outflow  radius  of  1.1  kpc  (see  Appendix  A  for  details).  We 
note  that  this  is  probably  a  conservative  estimate  of  the  outflow 
extension,  since  the  maps  of  the  wings  in  Fig.  [2jc-d)  suggest  the 
presence  of  even  more  extended  structures,  although  at  a  low  sig¬ 
nificance  level  (3cr).  The  molecular  outflow  discovered  in  IRAS 
F10565+2448  appears  to  be  less  powerful  than  in  the  other  Sy- 
ULIRGs  Mrk  231  and  IRAS  F08572+3915:  indeed  it  has  lower 
velocity  and,  carrying  a  mass  of  molecular  gas  of  ~  2  x  108  Me, 
has  a  mass-loss  rate  of  “only”  300  M0  yr-1. 


IRAS  23365+3604 

Previous  IRAM  PdBI  observations  of  the  CO(l-O)  transition  in 
IRAS  23365+3604  suggest  the  presence  of  a  compact  rotating 
ring  or  disk  ([Downes  &  Solomon|1998|),  responsible  for  the  ob¬ 
served  narrow  line  profile  with  a  FWZI  ~  500  km  s-1.  Our 
new  observations  recover  a  total  CO  flux  of  Sco,tot  =  (48.90 
±  0.20)  Jy  km  s-1  (Table [2},  which  is  60%  larger  than  the  IRAM 


PdBI  flux  measured  by  Downes  &  Solomon  (1998]),  but  consis¬ 
tent  with  the  IRAM  30m  flux  reported  by  the  same  authors  (see 
also  [Papadopoulos  et  al.|  (|2012b|)) . 

In  Fig.  [3]  we  present  the  continuum- subtracted  CO(l-O) 
spectrum,  extracted  from  an  aperture  with  a  diameter  of  5  arc¬ 
sec  (a)  and  from  the  central  pixel  (b).  The  continuum  emis¬ 
sion  was  evaluated  in  the  velocity  ranges  ve(-3700,  -1200)  km 
s-1  and  ve  (1200,  4000)  km  s-1.  The  5  arcsec  aperture  spec¬ 
trum  shows  wings  of  the  CO(l-O)  emission  line  which  extend 
to  about  ~  600  km  s-1  from  the  systemic  velocity,  as  well  as 
some  low  significance  emission  at  higher  velocities.  The  CO(l- 
0)  wings  are  instead  less  prominent  in  the  spectrum  extracted 
from  the  central  pixel.  The  cleaned  maps  of  the  blue  and  red- 
shifted  CO(l-O)  wings  (Fig.  pjc-d)),  integrated  within  the  veloc¬ 
ity  ranges  (-600,  -300)  km  s^and  (300,  600)  km  s-1,  show  that 
these  are  detected  at  4<x  and  8<x  significance,  respectively.  By 
combining  their  emission,  we  reach  a  signal-to-noise  of  9  (the 
combined  map  and  uv  plot  are  shown  in  the  Appendix,  Fig.  |A.5|) 


The  position- velocity  diagram  along  the  rotation  major  axis 


(Fig.  BTe-f))  reveals  that  such  molecular  gas  at  |v|  >  300  km  s 


Article  number,  page  6  of  26 


C.  Cicone  et  al.:  Massive  Molecular  Outflows  and  Evidence  for  AGN  Feedback  from  CO  Observations 


10 


5  <5R.A°  (") 


Fig.  3:  Continuum-subtracted  IRAM-PdBI  spectra,  maps  and  position- velocity  diagram  of  the  CO(l-O)  emission  line  of  IRAS 
23365+3604.  Panels  (a)  and  (b)  show  the  spectra  extracted  from  a  circular  aperture  with  diameter  of  5  arcsec  and  from  the  central 
pixel,  respectively.  For  display  purposes,  the  spectra  are  re-binned  in  channels  of  54  km  s-1.  In  both  panels  the  narrow  core  of 
the  line  is  fitted  with  a  single  Gaussian  function  (red  dashed  line),  (c-d)  Cleaned  map  of  the  emission  integrated  in  the  blue  and 
red-shifted  CO(l-O)  wings.  Contours  correspond  to  l<x  (l<x  rms  level  of  the  two  maps  is  0.2  mJy  beam-1).  The  physical  scale  at  the 
redshift  of  the  source  is  1.232  kpc  arcsec-1.  (e)  Position- velocity  diagram  along  the  major  axis  of  rotation.  The  inset  (f)  highlights 
the  rotation  pattern  traced  by  the  narrow  component  of  the  CO(l-O)  line;  this  also  suggests  the  presence  of  a  counter-rotating  disk 
visible  at  large  radii  (i.e.  up  to  about  8  arcsec  from  the  nucleus,  which  correspond  to  ~10  kpc).  Contours  are  in  steps  of  lcr  (up 
to  lcr)  in  panel  (e),  and  5 cr  in  the  inset  (f);  negative  contours  are  in  steps  of  lcr  (magenta  dashed  lines).  Each  pixel  corresponds  to 
0.6  arcsec  x  54  km  s-1  in  diagram  (e),  and  to  0.3  arcsec  x  27  km  s-1  in  the  inset  (f).  Note  the  signature  of  CO(l-O)  emission  at 
velocities  higher  than  ±300  km  s-1,  which  does  not  follow  the  rotation  pattern  and  is  therefore  interpreted  as  associated  with  the 
molecular  outflow  detected  by  Herschel-PACS. 


does  not  follow  the  rotation  pattern  traced  by  the  central  core 
emission:  this  suggests  that  this  gas  is  tracing  the  same  molecu¬ 
lar  outflow  discovered  by  Herschel-PACS  ( Veilleux  et  al.|2013|). 

It  is  interesting  to  note  that  the  position- velocity  diagram  ex¬ 
hibits  the  signature  of  a  second  disk  or  ring,  counter-rotating 
with  respect  to  the  bulk  of  the  molecular  gas,  which  was  not 
detected  by|Downes  &  Solomon|(|1998]).  This  may  provide  fur¬ 
ther  support  to  the  hypothesis  that  IRAS  23365+3604  is  a  later 
merger. 

We  estimate  the  molecular  outflow  in  IRAS  23365+3604  to 
carry  a  mass  of  molecular  gas  of  ~  1.5  x  108  M©  and  to  have  an 
approximate  radius  of  1.2  kpc  (obtained  by  fitting  directly  the 
uv  data),  resulting  in  an  outflow  mass-loss  rate  of  170  M©  yr-1, 
significantly  smaller  than  in  the  other  ULIRGs. 


Mrk  273 

The  continuum- subtracted  CO(l-O)  spectra  of  Mrk  273,  ex¬ 
tracted  from  a  circular  aperture  of  5  arcsec  diameter  and  from 
the  central  pixel,  are  shown  in  Figure  |4|a-b).  The  continuum 


emission  was  estimated  in  the  velocity  ranges  ve(-3000,  -2000) 
km  s-1  and  ve  (2000,  3700)  km  s-1.  The  total  CO(l-O)  line  flux 
that  we  measure  within  +  2000  km  s-1  is  Sco,tot  =(90.60  + 
0.50)  Jy  km  s-1  (Table [2},  consistent  with  the  flux  recovered  by 
previous  single  dish  observations  ([Sanders  et  al.|1991[|Solomon| 


et  al. 1 1997] | Albrecht  et  al.|2007 ;  Papadopoulos  et  al.  2012b),  and 


about  15%  larger  than  the  value  obtained  by|Downes  &  Solomon 
(1998)  with  the  PdBI.  The  CO(l-O)  line  spectra  show  a  broad 
emission  core,  which  has  already  been  detected  by  previous  ob¬ 
servations,  mostly  tracing  molecular  gas  in  a  nuclear  rotating 
disk  (see  Appendix  for  details).  We  fit  this  spectral  core  com¬ 
ponent  with  the  sum  of  two  Gaussian  functions  (red  dashed  line 
in  Fig.  Ef  -b).  In  addition  to  this  component,  our  new  obser¬ 
vations,  obtained  with  a  much  larger  bandwidth  than  previous 
data,  reveal  extra  CO(l-O)  emission  at  high  blue  and  red-shifted 
velocities.  Such  broad  blue  and  red  CO  wings  are  detected  at  a 
significance  of  6 cr  and  9 cr,  respectively,  as  shown  by  their  inte¬ 
grated  maps  in  panels  (c)  and  (d).  The  rotation  in  the  central 
concentration  of  molecular  gas  is  not  clearly  resolved  by  our 
D-configuration  (i.e.  low  spatial  resolution)  observations  (Fig. 
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Fig.  4:  Continuum- subtracted  IRAM-PdBI  spectra,  maps  and  position-velocity  diagram  of  the  CO(l-O)  emission  line  of  Mrk  273. 
Panels  (a)  and  (b)  show  the  spectra  extracted  from  a  circular  aperture  with  diameter  of  5  arcsec  and  from  the  central  pixel,  respec¬ 
tively.  For  display  purposes,  the  spectra  are  re-binned  in  channels  of  53  km  s-1  (and  in  channels  of  160  km  s-1  in  the  insets).  In 
both  panels  (a)  and  (b)  the  narrow  core  of  the  line  is  fitted  with  two  Gaussian  functions  (red  dashed  lines,  further  comments  in  the 
text  and  in  the  Appendix).  Note  the  appearance,  in  both  spectra,  of  broad  CO  wings,  and  in  particular  with  velocities  higher  than 
500  km  s-1.  (c-d)  Cleaned  map  of  the  emission  integrated  in  the  blue  and  red-shifted  CO(l-O)  wings.  Contours  correspond  to  Icr 
(1  cr  rms  noise  of  the  two  maps  is  0.2  mJy  beam-1).  The  physical  scale  at  the  redshift  of  the  source  is  0.746  kpc  arcsec-1.  (e-f) 
Position-velocity  diagram  along  the  major  axis  of  rotation.  Note  that  our  D-configuration  observations  do  not  completely  resolve 
the  disk  rotation  in  this  source.  Contours  are  in  steps  of  lcr  (up  to  10cr)  in  panel  (e),  and  5cr  in  the  inset  (f);  negative  contours  are  in 
steps  of  lcr  (magenta  dashed  lines).  Each  pixel  corresponds  to  2.2  arcsec  x  120  km  s-1  in  diagram  (e),  and  to  1.1  arcsec  x  54  km 
s-1  in  the  inset  (f).  Gas  at  velocities  higher  than  400  km  s-1  does  not  follow  the  rotation  pattern  and  is  most  likely  associated  to  the 
molecular  outflow  detected  by  Herschel-PACS. 


Ef-  However,  PdBI  observations  of  the  CO(2-l)  transition  pre¬ 
sented  by  Downes  &  Solomon  1998  resolve  the  kinematics  of 
the  nuclear  disk  of  Mrk  273,  revealing  a  velocity  centroid  that 
changes  from  -300  to  200  km  s-1.  The  PV  diagram  in  panels 
(e-f)  shows  that  the  CO(l-O)  emission  at  |v|  >  400  km  s-1  does 
not  follow  the  central  rotation  pattern,  and  is  therefore  ascribable 


to  the  molecular  outflow  detected  by  Herschel-PACS  (Veilleux 
eTaTpPB).  We  note  that  even  if  restricting  the  velocity  range 


of  integration  to  velocities  |v|  >  500  km  s  \  the  broad  CO(l-O) 
wings  of  Mrk  273  are  still  both  detected  at  a  significance  of  5<x. 


ifestly  visible  in  the  spectra  in  Fig.  |4ja-b),  is  tracing  the  low- 
velocity  component  of  a  prominent  red- shifted  outflow.  Such 
outflow  is  mostly  extended  to  the  North  up  to  velocities  of  500 
km  s-1,  while  it  is  more  compact  at  larger  velocities. 

For  the  molecular  outflow  in  Mrk  273,  we  estimate  an  aver¬ 
age  radius  of  about  600  pc  (see  Appendix  for  details)  and  a  mass 
of  ~  1.7  x  108  M©,  yielding  a  mass-loss  rate  of  600  M©  yr-1. 


IRAS  F23060+0505 


In  principle,  the  need  for  two  Gaussians  to  reproduce  the  core 
of  the  CO(l-O)  emission  is  not  an  anomaly  (see  also  the  case  of 
Mrk  231,  Cicone  et  al.  [2012|,  but  it  can  reflect  the  kinematics  of 
the  rotating  gas  in  the  central  molecular  disk.  However,  we  can¬ 
not  rule  out  the  possibility  that  the  asymmetry  on  the  red  side  of 
the  line  core,  which  is  fitted  by  the  smaller  Gaussian  centred  at 
about  250  km  s-1,  is  also  a  signature  of  outflow.  Indeed  we  de¬ 
tect  in  our  channel  maps  (not  shown)  a  blob  elongated  towards 
the  North  at  velocities  v  >  150  km  s-1,  further  confirming  our 
hypothesis  that  the  red- shifted  peak  of  CO(l-O)  emission,  man¬ 


To  our  knowledge,  this  is  the  first  published  observation  of  the 
CO(l-O)  molecular  transition  in  the  powerful  type  2  QSO  and 
ULIRG  IRAS  F23060+0505.  We  measure  a  total  integrated  CO 
flux  of  S co, tot  =  (15.30  ±  0.20)  Jy  km  s-1.  The  continuum- 
subtracted  CO(l-O)  spectra,  extracted  from  an  aperture  with  di¬ 
ameter  of  10  arcsec  and  from  the  central  pixel,  are  shown  in  pan¬ 
els  (a-b)  of  Figure  [5]  The  spectrum  extracted  from  a  10  arcsec 
aperture  shows  some  indication  of  broad  CO  emission,  which 
extends  up  to  at  least  v  ~  1000  km  s-1  in  the  red- shifted  side. 
Such  evidence  of  high  velocity  emission  is  much  more  marginal 
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Fig.  5:  Continuum-subtracted  IRAM-PdBI  spectra,  maps  and  position-velocity  diagram  of  the  CO(l-O)  emission  line  in  IRAS 
F23060+0505.  Panels  (a)  and  (b)  show  the  spectra  extracted  from  a  circular  aperture  with  diameter  of  10  arcsec  and  from  the 
central  pixel,  respectively.  For  display  purposes,  the  spectra  are  re-binned  in  channels  of  60  km  s-1.  (c-d)  Cleaned  maps  of  the 
emission  integrated  in  the  blue  and  red- shifted  CO(l-O)  wings.  Contours  correspond  to  l<x  (l<x  rms  noise  level  is  0.4  mJy  beam-1 
in  the  blue  wing  and  0.2  mJy  beam-1  in  the  red  wing).  The  physical  scale  at  the  redshift  of  the  source  is  2.934  kpc  arcsec-1.  (e-f) 
Position-velocity  diagram  along  the  major  axis  of  rotation.  Panel  (f)  highlights  the  rotational  kinematics  traced  by  the  narrow  core 
of  the  CO(l-O)  emission.  Contours  are  in  steps  of  ler  (up  to  10<x)  in  panel  (e)  and  of  5<x  in  the  inset  (f).  Negative  contours  are  in 
steps  of  1  cr  (magenta  dashed  lines).  Each  pixel  corresponds  to  0.7  arcsec  x  60  km  s-1  in  both  diagrams  (e)  and  (f).  Note  that  the 
CO(l-O)  wings  at  |v|  >  300  km  s-1  deviate  from  the  rotation  pattern,  and  are  likely  tracing  an  outflow  of  molecular  gas,  similarly  to 
the  one  discovered  in  the  other  sources. 


in  the  central  pixel  spectrum.  An  interesting  emission  feature  at 
red-shifted  velocities  of  v  e  (200, 400)  km  s-1  is  clearly  distin¬ 
guishable  in  both  spectra. 

The  cleaned  maps  of  the  broad  CO  wings,  integrated  within 
velocities  of  -500  <  v  <  -300  km  s-1  and  300  <  v  <  1100 
km  s-1,  are  presented  in  panels  (c-d)  of  Fig.  |5J  these  show 
that  the  blue  and  red  wings  are  detected  at  9cr  and  6cr  signifi¬ 
cance.  Moreover,  the  position- velocity  diagram  in  Fig.  [5je-f) 
reveals  that  such  emission  at  |v|  >  300  km  s-1  is  clearly  deviat¬ 
ing  from  the  central  rotational  pattern.  This  is  particularly  evi¬ 
dent  for  the  red-shifted  CO  emission:  the  PV  diagram  confirms 
that  the  feature  appearing  in  the  spectra  at  v  e  (200, 400)  km  s-1 
does  not  trace  rotating  gas,  but  it  is  likely  tracing  either  an  out¬ 
flow  or  an  inflow  of  molecular  gas.  The  blue- shifted  gas  at 
v  <  -300  km  s-1  is  mostly  elongated  to  the  north  of  the  galaxy, 
as  shown  by  the  map  in  panel  (c)  of  Fig.  [5]  We  also  note  the  de¬ 
tection,  in  the  PV  diagram  (as  well  as  in  the  10  arcsec  aperture 
spectrum),  of  very  high  velocity  (v  >  500  km  s-1)  red-shifted 
CO  emission,  although  at  a  low  significance  level  (i.e.  S/N<5). 
As  already  mentioned  in  Section  3.1,  according  to  our  criterion, 
IRAS  F23060+0505  does  not  fully  qualify  as  a  reliable  “out¬ 
flow”  detection,  since  high  velocity  gas  is  only  detected  at  low 


significance  and  we  lack  independent  confirmation  from  Her- 
schel  data  of  the  presence  of  a  molecular  outflow  in  this  object. 
Additional  observations  are  certainly  required  to  understand  the 
nature  of  the  high  velocity  CO  emission  in  this  powerful  QSO. 

We  use  the  flux  and  the  spatial  extent  of  the  broad  CO  wings 
(estimated  to  be  as  large  as  4  kpc  in  radius,  see  Appendix  for 
further  details)  to  estimate  upper  limits  on  the  outflow  mass 
and  mass-loss  rate  in  this  source,  which  may  be  as  large  as 
~  4  x  109  M©  and  1500  M©  yr-1,  respectively. 

Mrk  876  (PG  1613+658) 

Previous  interferometric  and  single  dish  observations  of  the  un¬ 
obscured  QSO  and  LIRG  Mrk  876  detected  the  CO(l-O)  emis¬ 
sion  at  high  signal-to-noise  ([Evans  et  al.||2001|  |2006),  and  es¬ 
timated  a  total  CO  integrated  flux  that  is  completely  consistent 
with  our  PdBI  measure  of  Sco,tot  =  (8.5  ±  1.5)  Jykms-1.  The 
continuum- subtracted  PdBI  CO(l-O)  spectrum,  extracted  from 
an  aperture  with  diameter  of  10  arcsec  and  from  the  central  pixel, 
is  presented  in  Fig.  |6ja-b).  The  signal-to-noise  of  our  spectrum 
is  low  due  to  the  short  integration  (6  hours  of  on  source  time  with 
the  PdBI,  see  Table  [I]).  Both  the  aperture  and  the  central  pixel 
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Fig.  6:  Continuum-subtracted  IRAM-PdBI  spectra  and  position-velocity  diagram  of  the  CO(l-O)  emission  line  of  Mrk  876.  Panels 
(a)  and  (b)  show  the  spectra  extracted  from  a  circular  aperture  with  diameter  of  10  arcsec  and  from  the  central  pixel,  respectively. 
For  display  purposes,  the  spectra  are  re-binned  in  channels  of  57  km  s-1.  No  clear  evidence  for  broad  wings  is  observed,  (c) 
Position-velocity  diagram  along  the  major  axis  of  the  CO(l-O)  rotation.  Contours  are  in  steps  of  lcr;  negative  contours  are  in  steps 
of  lcr  (magenta  dashed  lines).  Each  pixel  in  diagram  (c)  corresponds  to  1.3  arcsec  x  29  km  s-1.  The  physical  scale  at  the  redshift 
of  the  source  is  2.295  kpc  arcsec-1. 
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Fig.  7:  Continuum- subtracted  IRAM-PdBI  spectra  and  position- velocity  diagram  of  the  CO(l-O)  emission  line  of  I  Zw  1.  Panels  (a) 
and  (b)  show  the  spectra  extracted  from  a  circular  aperture  with  diameter  of  10  arcsec  and  from  the  central  pixel,  respectively.  For 
display  purposes,  the  spectra  are  re-binned  in  channels  of  54  km  s-1.  No  evidence  for  broad  wings  is  observed,  (c)  Position- velocity 
diagram  along  the  major  axis  of  rotation.  Contours  correspond  to:  lcr,  3cr,  6<x,  12cr,  18cr,  30cr,  36<x,  39<x,  42cr,  48cr;  negative 
contours  are  in  steps  of  lcr  (magenta  dashed  lines).  Each  pixel  in  diagram  (c)  corresponds  to  1.2  arcsec  x  27  km  s-1.  The  physical 
scale  at  the  redshift  of  the  source  is  1.175  kpc  arcsec-1. 
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spectra  do  not  show  any  clear  evidence  of  broad  CO  wings.  We 
marginally  detect  some  high  velocity  blue  and  red-shifted  CO 
emission  in  the  maps,  obtained  by  averaging  the  frequency  chan¬ 
nels  corresponding  to  velocities  of  (-500,  -300)  km  s-1  and  (400, 
1700)  km  s-1  (not  shown).  The  significance  of  this  detection 
is  however  low.  Moreover,  in  the  position-velocity  diagram  (in 
Fig.  [6]c)  there  is  no  apparent  trace  of  the  broad  CO  wings  that  we 
detect  in  the  maps.  For  this  reason,  and  because  Herschel-PACS 
observations  of  the  OH  transition  in  this  source  did  not  find  any 
evidence  for  a  blue-shifted  OH  absorption  (which  would  be  in¬ 
dicative  of  a  molecular  outflows),  Mrk  876  does  not  satisfy  our 
criterion  for  claiming  the  detection  of  a  massive  molecular  out¬ 
flow. 

By  combining  the  blue  and  red-shifted  CO(l-O)  emission, 
we  obtain  enough  S/N  (i.e.  S/N  =10,  Fig.  |A.11|)  to  evalu¬ 
ate  the  physical  extent  of  the  broad  wings  through  the  analysis 
of  the  uv  data  (see  Appendix  for  details).  Similarly  to  IRAS 
F23060+0505,  we  use  the  flux  and  radius  of  3.6  kpc  of  these  pu¬ 
tative  CO  wings  to  estimate  the  upper  limits  of  ~  3  x  109  Me  on 
the  mass  of  outflowing  gas  and  of  1800  M0  yr-1  on  the  outflow 
rate  in  Mrk  876.  We  stress  that  additional  interferometric  obser¬ 
vations  are  required  to  investigate  this  source  in  details,  and  to 
draw  more  definite  conclusions  about  the  presence  and  the  prop¬ 
erties  of  a  possible  molecular  outflow. 

Finally,  we  note  that  the  position- velocity  diagram  in  panel 
(c)  of  Fig.  [6]  shows  evidence  for  a  double  peak  (also  clearly 
discernible  in  the  spectrum  in  panel  (a)),  which  may  be  tracing  a 
ring-like  distribution  of  molecular  gas. 


I  Zw  1  (PG  0050+124) 


No  broad  CO(l-O)  emission  is  detected  in  the  spectra  and  maps 
of  the  LIRG  and  luminous  QSO  host  galaxy  I  Zw  1.  A  com¬ 
parison  between  our  PdBI  flux,  Sco,tot  =  (23.5  ±  1.0)  Jy  km 
s-1,  and  previous  single-dish  observations  of  this  galaxy  (Evans 


et  al.[||[2006t|Papadopoulos  et~al.  |2012b[),  shows  that  we  may  be 


missing  over  20%  of  more  extended  emission,  and  therefore  the 
presence  of  a  faint  and  broad  CO(l-O)  component,  cannot  be 
ruled  out  with  our  data.  We  note  that  this  object  does  not  show 
OH  absorption  either:  in  contrast,  OH  is  detected  in  pure  emis¬ 
sion,  red-shifted  by  ~  200  km  s-1  ( Veilleux  et  al.  2013).  We 
present  in  Fig.  |7Ja-b)  the  continuum- subtracted  CO(l-O)  spec¬ 
trum  of  I  Zw  1,  which  is  overall  narrow,  with  a  FWZI  of  500 
km  s-1.  The  position- velocity  diagram  obtained  along  the  major 
axis  of  the  CO(l-O)  rotation  (panel  (c)  of  Fig.  [7])  shows  a  clearly 
double-peaked  emission,  suggesting  that  I  Zw  1  may  host  a  ring¬ 
like  distribution  of  molecular  rotating  gas  in  its  nucleus. 


For  this  source,  we  estimate  upper  limits  for  the  flux  in  the 
broad  CO(l-O)  component  as  follows.  We  average  the  frequency 
channels  corresponding  to  the  velocity  ranges  of  ve(-650,  -250) 
km  s-1  and  ve(500,  750)  km  s-1,  and  produce  the  cleaned  maps, 
which  however  do  not  reveal  any  emission  above  the  noise  level. 
Our  upper  limits  for  the  fluxes  in  the  blue  and  red  CO  wings 
correspond  to  to  the  rms  noise  of  the  two  maps  (i.e.  0.13  mJy 
beam-1  and  0.19  mJy  beam-1)  multiplied  by  a  factor  of  3,  and 
integrated  over  400  and  250  km  s-1,  respectively.  We  use  these 
upper  limit  fluxes  and  an  hypothetical  CO  wings  radius  of  500 
pc  (i.e.  roughly  comparable  to  the  extent  of  the  narrow  core 
emission  in  I  Zw  1,  see  Appendix  for  more  details),  to  estimate 
upper  limits  to  the  outflowing  molecular  mass  and  outflow  rate 
in  this  source. 


4.  Extended  Sample 


4.1.  Description 

In  the  following  analysis  we  combine  our  7  sources  observed 
with  the  PdBI  with  a  heterogeneous  sample  of  12  local  galaxies 
in  which  outflows  of  molecular  gas  have  been  constrained  by 
mapping  the  two  lowest- J  CO  rotational  transitions  (i.e.  the  7=1- 
0  and  7=2-1).  The  complete  list  of  sources  is  given  in  Table  [3] 
In  the  following,  we  briefly  discuss  the  properties  of  the  sample 
of  sources  drawn  from  the  literature. 

The  Syl-ULIRG  Mrk  231  can  be  naturally  incorporated  in 
our  analysis,  since  its  powerful  molecular  outflow  was  studied 
using  the  same  methods  that  we  employ  for  the  new  sources  pre¬ 
sented  in  this  paper.  We  adopt  the  CO  integrated  fluxes  in  the 
blue  and  red-shifted  wings  and  the  outflow  size  (estimated  from 
the  7=1-0  transition)  that  are  listed  in|Cicone  et  al.|(|20T2]),  where 
the  PdBI  data  exploited  by|Feruglio  et  alA2010|  were  combined 
with  more  recent  observations. 

More  careful  consideration  is  needed  for  the  other  cases. 
|Alatalo  et  al.|(|201 1)  discovered  a  molecular  outflow  in  the  early- 
type  LINER-host  galaxy  NGC  1266,  using  a  combination  of  sin¬ 
gle  dish  (IRAM  30m)  and  interferometric  (SMA,  CARMA)  ob¬ 
servations  of  the  CO  7=1-0,  7=2-1  and  7=3-2  transitions.  To  cal¬ 
culate  the  outflow  mass,  we  exploit  the  flux  of  the  CO(l-O)  broad 
component  provided  by  these  authors.  By  assuming  the  bipo¬ 
lar  outflow  geometry  proposed  by|Alatalo  et  al.  ( 201 1|),  the  out¬ 
flow  size  (radius)  is  about  0.450  kpc,  which  refers  to  the  average 
width  of  the  blue  and  red-shifted  CO(2-l)  emission  as  measured 
with  the  SMA.  We  note  that,  by  using  the  7=2-1  transition,  we 


may  be  slightly  underestimating  the  outflow  size  (Cicone  et  al. 

|20l2|. 

The  molecular  outflow  in  the  starburst  galaxy  M82  has  been 
investigated  by  [Walter  et  al.|([2002|)  using  detailed  CO(l-O)  maps 
obtained  with  the  OVRO  interferometer.  We  adopt  the  CO  inte- 
grated  flux  of  Sco  =  7240  Jy  km  s-1,  estimated  by  Walter  et  al. 


(|2002|)  for  the  outflow  north  and  south  of  the  molecular  disk. 
The  wind  geometry  is  almost  spherical,  with  a  radius  of  about 
0.8  kpc. 

The  lenticular  galaxy  NGC  1377  and  the  LINER  NGC  1266 
share  some  similarities.  SMA  CO(2-l)  observations  of  NGC 
1377  revealed  the  presence  of  a  molecular  outflow,  traced  by 
broad  wings  of  the  emission  line  profile  (| Aalto  et  al.|2012b|).  In 
this  case,  both  the  flux  of  the  broad  component,  determined  by 
fitting  the  spectrum  with  two  Gaussians  (narrow  plus  broad),  and 
the  outflow  radius  of  0.2  kpc  (under  the  hypothesis  of  a  biconical 
outflow),  were  estimated  using  the  CO(2-l)  transition.  The  same 
caveats  as  for  NGC  1266  are  valid,  with  the  addition  that,  in  this 
object,  the  outflow  mass  is  also  obtained  using  the  CO(2-l)  line 
instead  of  the  CO(l-O)  line.  To  estimate  the  total  CO(l-O)  flux  in 
the  outflow  we  assume  thermalised  optically  thick  CO  emission, 
i.e.  L^O(10)  =  L^0(21)  (  Solomon  &  Vanden  Bout|2005 ). 

The  complex  morphology  of  the  CO(l-O)  emission  in  the 
AGN-dominated  ULIRG  and  merger  NGC  6240  has  been  in¬ 
vestigated  in  detail  using  the  compact  and  the  most  extended 
configurations  of  the  PdBI  ([Feruglio  et  al.[|2013a|b|).  |Feruglio[ 
|et  al.  p013b|)  clearly  show  that  the  blue-shifted  CO(l-O)  wing 
traces  molecular  gas  outflowing  from  the  southern  nucleus.  Con¬ 
versely,  most  of  the  red-shifted  emission  arises  from  the  region 
between  the  two  nuclei  and,  although  its  origin  is  strongly  de¬ 
bated,  it  is  believed  to  be  closely  related  to  the  merging  pro¬ 
cess.  Given  the  very  high  velocities  observed  in  this  red-shifted 
component,  it  is  tempting  to  presume  that,  among  the  other 
mechanisms  proposed  and  discussed  in  Feruglio  et  al.|  (|2013b|) 
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Object 

Total  Line 

Blue  Wing 

Red  Wing 

Wings  Size 

Sco 

Vel.  Range 

Sco 

Vel.  Range 

Sco 

FWHM 

(Jy  kms"1) 

(km ! 

S’1) 

(Jy  km  s-1) 

(km  s"1) 

(Jy  km  s-1) 

(kpc) 

IRAS  F08572+3915 

(10.80  +  0.70) 

(-1200, 

-400) 

(1.69  +  0.42) 

(500,  1100) 

(1.22  +  0.36) 

(1.64  +  0.34) 

IRAS  F10565+2448 

(108.0+  1.0) 

(-600, 

-300) 

(1.72  +  0.54) 

(300,  600) 

(1.50  +  0.24) 

(2.19  +  0.22) 

IRAS  23365+3604 

(48.90  +  0.20) 

(-600, 

-300) 

(0.40  +  0.12) 

(300,  600) 

(0.46  +  0.12) 

(2.45  +  0.70) 

Mrk  273 

(90.60  +  0.50) 

(-800, 

-400) 

(1.04  +  0.50) 

(400,  900) 

(2.07  +  0.52) 

(1.10  +  0.50) 

IRAS  F23060+0505 

(15.30  +  0.20) 

(-500, 

-300) 

(0.75  +  0.18) 

(300,  1100) 

(1.51  +  0.40) 

(8.1  ±2.9) 

Mrk  876 

(8.5  +  1.5) 

(-500, 

-300) 

(0.65  +  0.63) 

(400,  1700) 

(3.1  +  1.5) 

(7.1  +  1.4) 

IZw  1 

(23.5  +  1.0) 

(-650, 

-250) 

<0.16 

(500,  750) 

<0.14 

Non  detection 

and  references  therein,  there  may  be  an  outflow  contribution, 
but  the  current  data  do  not  allow  us  to  quantify  it.  We  there¬ 
fore  follow  Feruglio  et  al.|(|2Q13b])  and  estimate  the  outflow  rate 
from  the  blue  wing  of  the  CO(l-O)  line,  in  the  velocity  range 
ve(-500,  -200)  km  s-1,  which  has  an  integrated  flux  of  17.8 
Jy  km  s-1  and  a  total  size  of  1.3  kpc.  We  note  that  in 


Fer- 


|uglio  et  ah]  (|2013b|)  a  conservative  CO-to-H2  conversion  factor 
of  0.5  M0(Kkms  1  pc2)-1  has  been  used,  while  here  we  adopt 
the  same  value  of  0.8  as  in  the  other  ULIRGs. 

|Sakamoto  et  al.|  ([2Q06b|)  discovered  a  wind  of  molecular  gas 
in  the  starburst  merger  NGC  3256,  by  using  SMA  observations 
of  the  CO(2-l)  emission  line.  They  integrate  the  blue  and  red- 
shifted  high-velocity  emission  in  the  CO(2-l)  spectrum,  within 
the  velocity  ranges  (-195,  -165)  km  s-1  and  (215,  425)  km  s_1. 
The  resulting  total  CO(2-l)  flux  in  the  broad  wings  is  28  Jy  km 
s-1,  and  we  use  this  flux  to  estimate  the  mass  of  the  molecular 
outflowing  gas.  We  also  adopt  the  same  biconical  geometry  of 
the  outflow  as|Sakamoto  et  al.|(|2006b|),  with  an  outflow  exten¬ 
sion  of  0.5  kpc. 

The  starburst  galaxy  NGC  3628  hosts  a  sub-kpc  scale  molec¬ 
ular  outflow  to  the  north  of  the  galactic  disk  traced  by  CO(l-O) 
emission  (Tsai  et  al.|2012|).  This  is  spatially  coincident  with  the 
northern  ejection  point  of  a  large  scale  (~10  kpc)  plasma  out¬ 
flow,  which  instead  extends  both  to  the  north  and  to  the  south 
of  the  disk,  as  evidenced  by  soft  X-ray  Chandra  observations. 
Tsai  et  al.  (2012)  infer  a  CO(l-O)  flux  of  82  Jy  km  s-1  for 


the  molecular  outflow,  by  using  detailed  interferometric  (NMA) 
maps.  They  suggest  that  this  gas  is  expanding  isotropically  to 
the  northern  side  of  the  galactic  disk,  with  an  average  velocity  of 
50  km  s-1 ,  up  to  a  radius  of  0.4  kpc. 

The  CO(2-l)  emission  in  the  prototypical  starburst  galaxy 
NGC  253,  imaged  with  the  SMA  by|Sakamoto  et  al.|  |2006a|), 
displays  two  molecular  gas  features  at  a  distance  of  0.2  and 
0.7  kpc  from  the  galactic  centre,  which  have  been  interpreted 
as  expanding  super-bubbles.  Sakamoto  et  al.|  ([2006a)  estimate 
their  molecular  gas  masses  and  diameters  of  -100  pc  from  the 
CO(2-l)  channel  maps,  and  assume  cylindrical  symmetry  and 
non-linear  expansion  (i.e.  R(t)  oc  t1/2)  to  model  the  outflow.  We 
adopt  the  same  assumptions  as  these  authors,  but  we  apply  a  cor¬ 
rection  to  the  mass  estimates  to  account  for  the  use  of  a  different 
CO-to-H2_conversion  factor  (see  notes  in  Table]?]).  Very  recently, 


Bolatto  et  al.  (2013)  have  studied  with  unprecedented  detail  the 


molecular  outflow  of  NGC  253,  using  new  ALMA  observations 
of  the  CO(l-O)  transition.  The  outflow  rate  inferred  from  these 
new  data  of  M0f  ~  3  -  9  M©  yr_1 ,  is  in  excellent  agreement  with 
our  estimate  based  on  the  work  by|Sakamoto  et  al.|(|2006a|  (Ta- 
ble|§. 

The  integrated  CO(l-O)  PdBI  map  of  the  LINER  and  star- 
burst  galaxy  NGC  6764  shows  evidence  for  a  molecular  outflow 
extending  to  the  north  of  the  galactic  disk  (Leon  et  al.|2007|.  To 


estimate  the  outflow  rate,  we  use  the  CO  flux  of  1.4  Jy  km  s 


provided  by  [Leon  et  al.|(|2007|)  and  we  assume  that  this  is  emitted 
by  gas  uniformly  distributed  in  a  conical  volume  with  radius  of 
0.6  kpc,  outflowing  at  a  velocity  of  170  km  s-1. 


Multi-transitional  CO  interferometric  observations  have 
shown  the  presence  of  a  massive  molecular  outflow  in  the  Sy2- 
LIRG  NGC  1068,  revealed  by  broad  wings  of  the  molecular 
emission  lines  ([Krips  et  al.|2011|).  The  blue  and  the  red-shifted 
high  velocity  CO  components  are  co-spatial,  and  they  reach  ve- 
locity_shifts  of  up  to  250  km  s-1,  with  respect  to  the  systemic. 
|Krips  et  aL](2011|  suggest  that  up  to  30%  of  the  observed  PdBI 
CO(2-l)  emission  is  ascribable  to  gas  blowing  radially  outward 
from  the  galaxy  disk,  and  we  use  this  estimate  to  evaluate  the 
mass  of  the  outflowing  gas.  We  model  the  outflow  using  a  bicon¬ 
ical  geometry  with  radius  of  1  arcsec  (0.1  kpc)  and  an  average 
velocity  of  150  km  s-1. 

In  the  early  type  galaxy  IC  5063,  which  hosts  a  radio-loud 
Sy2  nucleus,  |Morganti  et  al.|(|2013|)  detect  with  APEX  a  promi¬ 
nent  blue-shifted  wing  of  the  CO(2-l)  line,  tracing  gas  with  ve¬ 
locities  exceeding  those  of  the  rotating  gas  in  the  disk.  Inter¬ 
estingly,  such  velocities  are  consistent  with  those  found  in  blue- 
shifted  HI  absorption,  also  detected  in  IC  5063,  tracing  an  out¬ 
flow  along  the  line  of  sight  of  the  radio  jet.  Since  the  APEX 
observations  do  not  provide  any  information  about  the  spatial 
extension  of  the  CO(2-l)  emission,  the  outflow  rate  estimation 
by|Morganti  et  al.  (  2013)  relies  on  the  (likely  but  not  yet  proven) 
assumption  that  the  molecular  and  neutral  winds  are  co-spatial. 
Therefore,  following  these  authors,  we  model  the  blue-shifted 
CO(2-l)  wing  as  a  conical  outflow  with  radius  of  0.5  kpc  and 
velocity  of  400  km  s-1.  Similarly  to  the  other  cases,  we  assume 


I '  -  I ' 

^CO(l-Q)  ~  lO(2-1)’ 

|Tsai  et  alT|  ([2009])  discovered  various  diffuse  molecular  fea¬ 
tures  extended  above  and  below  the  molecular  disk  of  the  nearby 
star-bursting  galaxy  NGC  2146,  by  using  very  deep  NMA  CO(l- 
0)  observations.  One  of  these  structures,  which  can  be  clearly 
seen  in  the  PV  diagram  and  in  the  integrated  intensity  and  ve¬ 
locity  maps,  traces  a  molecular  outflow  with  velocities  of  up  to 
200  km  s-1,  extending  to  a  distance  of  about  25  arcsec  from 
the  galactic  centre  (corresponding  to  1.6  kpc  with  the  cosmol¬ 
ogy  adopted  in  this  paper).  The  CO(l-O)  integrated  flux  in  this 
component  is  150  Jy  km  s-1. 


4.2.  Outflow  Properties 

As  for  the  seven  sources  of  our  PdBI  sample,  to  estimate  the  to¬ 
tal  mass  of  outflowing  molecular  gas  in  the  12  sources  belonging 
to  the  literature  sample,  we  adopt  a  unique  CO-to-H2  conver¬ 
sion  factor,  i.e.  the  value  commonly  used  for  ULIRG-like  gas 
conditions  (aco(i-o)  =  0-8  M©  (Kkm  s-1pc2)-1,  see  also  expla¬ 
nation  in  Section  3.2).  However,  for  the  non-ULIRG  galaxies 
belonging  to  this  sample,  we  also  report,  in  Table  [4]  the  out- 
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flow  mass,  mass-loss  rates  and  energetics  as  obtained  by  using 
the  same  orco(i-o)  adopted  by  the  authors  of  the  respective  pa¬ 
pers  for  estimating  the  total  amount  of  molecular  gas  in  that  spe¬ 
cific  host  galaxy  (see  column  6  of  Tableland  notes  in  Table [4]). 
When  only  the  CO(2-l)  transition  is  available  for  estimating  the 
outflow  mass,  as  in  the  case  of  NGC  1377  and  NGC  3256,  we 
assume  that  the  CO  emission  is  thermalised  and  optically  thick, 
hence  the  line  luminosity  L^n  is  independent  of  the  transition  J 
((Solomon  &  Vanden  Bout|2005]). 

In  the  12  sources  from  the  literature,  we  derive  the  molecular 
outflow  properties  by  assuming  the  same  geometry  as  in  the  cor¬ 
responding  papers.  These  works  usually  provide  the  dynamical 
time-scale  r^yn  of  the  wind,  and  calculate  the  outflow  mass-loss 
rate  using  the  relation:  M0f  -  M0f /Tdyn •  As  already  explained 
in  Sect.  3.2,  this  method  is  appropriate  in  a  shell-like  explosive 
scenario,  where  one  or  more  clouds  are  clearly  resolved  and  their 
distance  from  the  centre  can  be  measured.  However,  if  the  out¬ 
flowing  gas  appears  to  be  uniformly  distributed  within  a  spher¬ 
ical  or  multi-conical  volume,  our  prescription  is  more  correct; 
hence,  in  the  latter  case,  we  re-calculate  the  outflow  rate  with 
our  method. 


5.  Ancillary  Information 

In  this  section  we  discuss  how  the  ancillary  information  listed 
in  Table  3  was  derived.  The  star  formation  rates  in  column  4 
are  inferred  from  the  total  IR  (8-1000  /im)  luminosity  by  assum¬ 
ing  the  relation  in  Kennicutt  (1998  )  and  a  Chabrier  IMF  (see 
also  Sturm  et  al.  (2011 )):  SFR  =  (1  -  a)  •  10_10LIR,  using  the 


1041  SFR(M0  yr-1)  (see  also 


Maiolino  et  al. 


factor  a  =  Lagn/LboI  reported  in  column  9.  From  these  SFRs, 
we  calculate  the  kinetic  power  injected  by  supernovae  (column 
5),  by  using  the  relationship  provided  by  [Veilleux  et  aL  (2005): 
Pkin,SF(ergs_1)  =  7  . . 

|20l2l>. 

The  total  mass  of  molecular  gas  in  column  7  is  corrected  for 
He  and  is  derived  from  the  CO  integrated  flux,  using,  for  each 
source,  the  corresponding  CO-to-H2  conversion  factor  in  column 
6.  The  value  of  the  CO-to-H2  conversion  factor  in  different  envi¬ 
ronments  and  its  dependance  on  the  galaxy  properties  is  highly 
debated,  and  it  may  constitute  a  major  source  of  uncertainty  in 
our  estimates  of  the  molecular  gas  masses.  However,  such  dis¬ 
cussion  goes  beyond  the  scope  of  this  paper.  Here  we  simply  as¬ 
sume  arcoa-o) 
in  M82  and  M82 


=  0.8  in  ULIRGs  (|Downes  &  Solomon|1998|,  1.2 
[82-like  galaxies  (|Weifi  et  al.  2001),  4.36  (rounded 


to  4.4  in  Table  3)  in  the  other  galaxies  (Milky  Way  value,  Genzel 
[et  al.|2012|). 

Obtaining  a  reliable  estimate  of  the  AGN  luminosity  and  of 
the  AGN  contribution  to  the  total  bolometric  luminosity  of  a 
galaxy  is  in  general  a  difficult  task,  but  it  is  even  more  com¬ 
plicated  in  the  case  of  the  heavily  obscured  objects  included  in 
our  sample.  For  the  AGN-host  (U)LIRGs,  we  mostly  rely  on 
the  analysis  of  their  infrared  properties  as  revealed  by  Spitzer- 
IRS  observations;  in  particular,  for  IRAS  F08572+3915,  IRAS 
23365+3604,  Mrk  273,  IRAS  F23060+0505  and  Mrk  231,  we 
adopt  the  Lagn/LboI  ratios  inferred  using  a  spectral  decom¬ 
position  method,  which  relies  on  the  use  of  AGN  and  star- 
burst  templates,  in  the  rest- frame  5-8  /im  wavelength  range, 
where  the  starbursts  and  AGNs  have  quite  distinct  spectral 
properties  ([Nardini  et  al.||2009l  |2010|).  In  the  case  of  IRAS 
F10565+2448,  we  assume  the  average  AGN  bolometric  frac¬ 
tion  measured  by  Veilleux  et  al.|  ([2009),  using  a  combination 
of  continuum,  emission  and  absorption  line  diagnostics  based 
on  MIR  and  FIR  Spitzer  observations  of  this  galaxy.  We  then 
evaluate  for  these  ULIRGs  the  AGN  bolometric  luminosity  by 


assuming  LBoi  =  1.15  Lir  ([Veilleux  et  al.||2009|).  In  the  case  of 
I  Zw  1,  Mrk  876,  NGC  6240,  NGC  6764  and  IC  5063,  we  esti¬ 
mate  LAgn  from  their  absorption-corrected  hard  X-ray  luminos¬ 
ity  Lx (2  -  lOkeV)  using  the  relation  in  [Marconi  et  al.|(|2004). 
For  the  other  two  LINERs,  NGC  1266  and  NGC  1377,  we  calcu¬ 
late  LAgn  rising  the  [OIII]/15007  luminosity  (reported  by|Mous-] 
takas  &  Kennicutt  2006  )  as  a  tracer  of  AGN  activity,  following 
Heckman  et  al.|  ([2004 ).  For  the  extremely  obscured  Seyfert  2 


NGC  1068,  we  adopt  directly  the  AGN  luminosity  calculated  by 
[Prieto  et  aL|([2010|).  For  the  starburst  galaxies  M82,  NGC  3256, 
NGC  3628,  NGC  253  and  NGC  2146  ,  we  assume  Lb0i  -  Lir 
and  set  as  an  upper  limit  for  the  LAgn  the  value  estimated  from 
the  unabsorbed  hard  X-ray  luminosity  Lx(2  -  10  keV),  using  the 
relation  provided  by  Marconi  et  al.  (|2004|). 


6.  Discussion 

Recent  works  based  on  the  analysis  of  Herschel-PACS  FIR  ob¬ 
servations  of  local  ULIRGs  and  AGN  hosts  have  shown  that 
molecular  outflows  are  rather  common  in  this  class  of  objects 
([Veilleux  et  al.|2013}|Spoon  et  al.|2013).  In  our  study,  thanks  to 
the  additional  information  provided  by  the  interferometric  CO 
observations,  we  find  that  molecular  outflows  can  be  very  mas¬ 
sive  and  energetic,  with  outflow  rates  of  several  100  M©  yr-1, 
and  that  the  outflowing  gas  extends  on  kiloparsec  scales,  hence 
affecting  the  host  galaxy. 

There  is  a  good  correspondence  between  the  detection  of 
molecular  outflows  through  the  OH  P-Cygni  profile  and  the  de¬ 
tection  of  the  same  wind  through  the  CO  wings.  In  particular, 
we  note  that  in  the  case  of  IRAS  08572+3915  (as  well  as  for 
Mrk  231  and  for  the  starburst  galaxy  NGC  253),  the  outflow 
mass-loss  rate  estimate  obtained  by|Sturm  et  al.||201l|  by  mod¬ 
elling  the  OH  P-Cygni  profile  is  in  good  agreement  with  our 
results.  For  the  other  three  sources  in  which  OH  P-Cygni  pro¬ 
files  were  detected  by  Herschel  and  which  we  followed  up  with 
the  PdBI  in  their  CO(l-O)  transition  (i.e.  IRAS  10565+2448, 
IRAS  23365+3604  and  Mrk  273),  we  find  a  good  correspon¬ 
dence  between  the  outflow  velocities  traced  by  the  CO  broad 
wings  and  the  velocity  of  the  OH  blue-shifted  absorption  re¬ 
ported  by  [Veilleux  et  al.|([2013|.  Our  PdBI  observations  of  Mrk 
876  and  I  Zw  1  only  allow  us  to  set  upper  limits  on  the  molec¬ 
ular  outflow  properties  in  these  sources.  These  two  galaxies 
share  some  similarities:  both  are  dominated  by  a  powerful  unob¬ 
scured  QSO,  and  both  exhibit  pure  OH  emission  profiles  in  the 
Herschel-PACS  data  presented  by  Veilleux  et  al.  |2013|,  with  no 
evidence  for  blue- shifted  absorption.  However,  the  OH  profile 
of  Mrk  876  is  the  broadest  among  the  sources  with  OH  detected 
exclusively  in  emission:  this  suggests  non-gravitational  motion 
of  molecular  gas.  The  poor  signal-to-noise  of  Herschel-PACS 
observations  of  Mrk  876  and  I  Zw  1  does  not  allow  us  to  draw 
any  conclusion  about  the  possibility  of  very  high-velocity  (i.e. 
|v|  >  1000  km  s"1)  molecular  outflows. 

In  the  following  we  discuss  the  relationships  between  the 
properties  of  the  molecular  outflows  and  the  physical  properties 
of  the  AGN  and  of  the  host  galaxies,  to  shed  light  on  the  ori¬ 
gin  and  nature  of  the  molecular  outflows,  and  their  role  in  the 
evolution  of  galaxies.  A  general  caveat  of  this  study  is  that  our 
sample  is  certainly  not  representative  of  the  local  population  of 
galaxies.  It  is  incomplete  and  subject  to  biases.  Indeed,  the  bulk 
of  the  SHINING  subsample  (4  targets  out  of  5)  has  been  selected 
to  show  evidence  of  molecular  outflows  according  to  the  OH  P- 
Cygni  profiles.  As  a  consequence,  only  IRAS  F23060+0505, 
Mrk  876  and  I  Zw  1  were  observed  with  the  PdBI  without  a 
previous  knowledge  of  the  presence  of  a  molecular  outflow.  In 
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Table  3:  Properties  of  the  Extended  Sample  of  Galaxies 


Object 

type 

zCO 

SFR 

log(Pkin,SF) 

tfCO(l-O) 

log(Mf[2TOT  )' 

log(LAGN) 

Lagn/LboI 

Refs. 

[Moyr-1] 

[erg  s-1] 

[M0/(Kkms-1  pc2)] 

[Me] 

[erg  s-1] 

(1) 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

(8) 

(9) 

(10) 

IRAS  F08572+3915 

Sy  2 

0.05821 

20 

43.15 

0.8 

9.18 

45.72 

0.860 

1,  b 

IRAS  F10565+2448 

Sy  2 

0.04311 

95 

43.82 

0.8 

9.90 

44.81 

0.170 

1,  c 

IRAS  23365+3604 

LINER 

0.06438 

137 

43.98 

0.8 

9.93 

44.67 

0.072 

1,  b 

Mrk  273 

Sy  2 

0.03777 

139 

43.99 

0.8 

9.70 

44.73 

0.080 

1,  a 

IRAS  F23060+0505 

Sy  2 

0.17300 

75 

43.72 

0.8 

10.39 

46.06 

0.780 

1,  b 

Mrk  876 

Sy  1 

0.12900 

6.5 

42.66 

0.8 

9.84 

45.84 

0.930 

1,  d 

IZw  1 

Sy  1 

0.06114 

36 

43.40 

0.8 

9.56 

45.37 

0.520 

1,  d 

Mrk  231 

Sy  1 

0.04217 

234 

44.21 

0.8 

9.73 

45.72 

0.340 

2,  b 

NGC  1266 

LINER 

0.00719 

1.6 

42.05 

4.4 

9.23 

43.31 

0.250 

3,  f 

M  82 

HII 

0.00068 

10 

42.85 

1.2 

8.64 

<41.54 

<0.00090 

4,  i 

NGC  1377 

LINER 

0.00578 

0.9 

41.80 

4.4 

8.44 

42.93 

0.200 

5,  f 

NGC  6240 

Sy  2 

0.02448 

16 

43.05 

0.8 

9.86 

45.38 

0.780 

6,  e 

NGC  3256 

HII 

0.00926 

36 

43.40 

1.2 

9.68 

<41.97 

<0.00070 

7,j 

NGC  3628 

HII 

0.00280 

1.8 

42.10 

1.2 

9.53 

<40.79 

<0.00090 

8,  k 

NGC  253 

HII 

0.00081 

3 

42.32 

1.2 

8.15 

<40.66 

<0.00040 

9,j 

NGC  6764 

LINER 

0.00807 

2.6 

42.26 

1.2 

8.90 

42.23 

0.017 

10,  g 

NGC  1068 

Sy  2 

0.00379 

18 

43.10 

0.8 

9.11 

43.94 

0.097 

11, h 

IC  5063 

Sy  2 

0.01100 

0.6 

41.62 

4.4 

8.85 

44.30 

0.900 

12,1 

NGC  2146 

HII 

0.00298 

12 

42.92 

1.2 

8.94 

<41.09 

<0.00030 

13,  m 

The  following  information  is  listed  for  each  galaxy  in  each  of  the  columns:  (1)  name;  (2)  optical  classification;  (3)  redshift  based  on  the  CO 
spectrum;  (4)  star  formation  rate;  (5)  kinetic  energy  associated  with  the  supernovae;  (6)  CO-to-H2  conversion  factor  adopted  to  infer  MH2  from 
the  CO(l-O)  luminosity;  (7)  total  mass  of  molecular  gas;  (8)  AGN  bolometric  luminosity;  (9)  ratio  between  the  bolometric  luminosity  of  the  AGN 
and  the  total  bolometric  luminosity  of  the  galaxy;  (10)  references  for  the  CO  fluxes  used  to  estimate  MH 


-  Alatalo  et  al.  (2011),  4  -  Walter  et  al.|(2002),  5  -|Aalto  et  al.  f2012b),  6  -jFeruglio  et  al.  (2013a[b),  7  -|Sakamoto  et  al.  (2006b),  8  -  Tsai  et  al. 

(2012 

),  9  -  Mauersberger  et  al.  ( 

1996),  10  -panders  &  Mirabel  (1985  ),  11  -  Maiolino  et  al.  (1997  ),  12  -  Wiklind  et  al.  (1995  ),  13  -  Tsai  et  al. 

(2009 

);  references  for  LAgn-  a  - 

Nardini  et  al.  (2009  ),  b^lNardini  et  al.||2010),  c  -  Veilleux  et  aL)j(i2009  ),  d  -  Piconcelli  et  al.  (2005),  e  -  Risaliti 

(|2003|),  k  -|Gonzalez-Martm  et  al.|(|2009]>7l  -|Koyama  et  al.|(|1992|),  m  -|della  Ceca  et  al.|(|1999|>,|Inui  et  al.|p005|>. 

Notes:  t  In  these  cases  of  NGC  1377  and  NGC  3256,  we  estimate  the  total  molecular  gas  mass  by  using  the  CO(2-l)  transition  (instead  of  the 
CO(l-O)),  by  assuming  that  the  CO  gas  is  thermalised  and  optically  thick,  i.e.  that  L^0(1  0)  =  L^0(2 


Table  4:  Outflow  Properties  of  the  Extended  Sample  of  Galaxies 


Object 

log(Mn2,oF)T 

Mh2,oft 

Rof 

vOF,avg 

vOF,max 

log(Tdep)1' 

log(Pkin,OF)T 

(Mh2,of  v)/ (LAgn  /  c)T 

Refs. 

[Mo] 

[M0  yr1] 

[kpc] 

[km  s-1] 

[km  s-1] 

[yr] 

[ergs-1] 

(1) 

(2) 

(3) 

(4) 

(5) 

(6) 

(V) 

(8) 

(9) 

(10) 

IRAS  F08572+3915 

8.61 

1210 

0.82 

800 

1200 

6.10 

44.39 

35 

1 

IRAS  F10565+2448 

8.37 

300 

1.10 

450 

600 

7.42 

43.28 

40 

1 

IRAS  23365+3604 

8.17 

170 

1.23 

450 

600 

7.70 

43.04 

31 

1 

Mrk  273 

8.24 

600 

0.55 

620 

900 

6.92 

43.86 

130 

1 

IRAS  F23060+0505 

<9.56 

<  1500 

<4.05 

(550) 

(1100) 

>7.21 

<44.16 

<14 

1 

Mrk  876 

<9.48 

<  1830 

<3.55 

(700) 

(1700) 

>6.57 

<  44.45 

<35 

1 

IZw  1 

<7.67 

<  140 

(0.50) 

(500) 

(750) 

>7.41 

<  43.04 

<6 

1 

Mrk  231 

8.47 

1050 

0.60 

700 

1000 

6.71 

44.21 

26 

2 

NGC  1266 

7.93  -  8.66 

33  -  180 

0.45 

177 

362 

7.71  -  6.97 

41.51-42.25 

54  -  300 

3 

M  82 

8.08  -  8.25 

12-18 

0.80 

100 

230 

7.56-7.38 

40.58  -  40.75 

>650  -  >980 

4 

NGC  1377 

7.29  -  8.03 

14-76 

0.20 

110 

140 

7.29  -  6.56 

40.73  -41.46 

34  -  190 

5 

NGC  6240 

8.61 

800 

0.65 

400 

500 

6.96 

43.61 

25 

6 

NGC  3256 

7.34-7.51 

11  -  16 

0.50 

250 

425 

8.64  -  8.48 

41.34-41.50 

>560 ->810 

7 

NGC  3628 

7.36  -  7.54 

4.5  -  6.7 

0.40 

50 

110 

8.88  -  8.70 

39.55  -  39.72 

>690 ->1030 

8 

NGC  253 

6.32  -  6.50 

4.2 -6.3 

0.20 

50 

100 

7.53-7.35 

39.52  -  39.70 

>870 ->1300 

9 

NGC  6764 

6.52  -  6.69 

3.1  -4.7 

0.60 

170 

280 

8.41  -  8.23 

40.45  -  40.63 

59-89 

10 

NGC  1068 

7.26 

84 

0.10 

150 

250 

7.19 

41.77 

27 

11 

IC  5063 

7.37-8.10 

23  -  127 

0.50 

300 

450 

7.48  -  6.75 

41.82-42.56 

7-36 

12 

NGC  2146 

7.68-7.86 

14-22 

1.55 

150 

200 

7.78-7.61 

41.00-41.18 

>3300  -  >4900 

13 

The  following  information  is  listed  for  each  galaxy  in  each  of  the  columns:  (1)  name;  (2)  mass  of  the  molecular  gas  in  the  outflow;  (3)  molecular 
outflow  mass-loss  rate;  (4)  radius  of  the  outflow  extension;  (5)  average  velocity  of  the  outflow;  (6)  maximum  velocity  of  the  outflow;  (7) 
depletion  time-scale  associated  with  the  outflow;  (8)  kinetic  power  of  the  outflow;  (9)  outflow  momentum  rate  relative  to  LAGn/c;  (10)  references 


for  the  molecular  outflow  measurements:  1  -  this  work,  2  -  Cicone  et  al.| ( 2012 

l,  3  -|Alatalo  et  al.  (2011),  4  -  Walter  et  al.  |2002  ),  5  -  Aalto  et  al. 

(2012b  ),  6  -  Feruglio  et  al.  (2013b  ),  7  -  Sakamoto  et  al.  (  2006b 

),  8  -|Tsai  et  al. 

(2012),  9  -  Sakamoto  et  al.  (2006a)  and  Bolatto  et  al.  (2013),  10  - 

[Leon  et  al.  (2007),  1 1  -|Krips  et  al.  (2011),  12  -  Morganti  et  al. 

p013p3-  Tsai  etal.  (2009  ) 

Notes:  T  For  non-ULIRG  galaxies  we  report  two  values  of  the  outflow  mass,  outflow  rate,  depletion  time,  kinetic  power  and  momentum  rate  of 
the  outflow  (columns  2,  3,  7,  8,  9),  obtained  by  adopting  in  one  case  the  conservative  CO-to-H2  conversion  factor  ofCo(i-o)=0.8,  which  was  found 
appropriate  for  Mrk231,  and  in  the  other  case  the  same  arCo(i-o)  (column  6  in  table [3] )  as  for  estimating  the  total  molecular  gas  mass  in  the  same 
host  galaxy  (column  7  in  table [3J .  The  two  values  of  aCo  yield  respectively  to  the  lower  and  the  higher  estimates  of  MH2,of,  MH2  and  Pkin,  and 
conversely  for  Tdep  (see  also  discussion  in  the  text). 
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C.  Cicone  et  al.:  Massive  Molecular  Outflows  and  Evidence  for  AGN  Feedback  from  CO  Observations 


IRAS  F23060+0505  we  detect  molecular  gas  at  velocities  larger 
than  300  km  s-1  whose  kinematics  deviates  from  the  central  ro¬ 
tational  pattern  (traced  by  the  narrow  core  of  the  CO  line):  how¬ 
ever,  the  lack  of  detection  of  high  velocity  CO  emission  at  high 
significance  prevents  us  from  drawing  firm  conclusions  about  the 
outflowing  origin  of  this  gas.  As  mentioned  before,  in  Mrk  876 
and  I  Zw  1  we  do  not  find  any  significant  evidence  for  the  pres¬ 
ence  of  a  massive  molecular  outflow,  although  in  Mrk  876  we 
tentatively  detect  some  high  velocity  CO  emission.  Summaris¬ 
ing,  none  of  our  three  “unbiased”  sources  constitutes  a  reliable 
molecular  outflow  detection.  On  the  other  hand,  the  sample 
collected  from  the  literature  only  contains  sources  for  which  a 
molecular  outflow  has  been  constrained:  this  unavoidably  biases 
the  sample  towards  galaxies  for  which  a  molecular  outflow  has 
been  detected,  since  in  most  cases  in  which  an  outflow  has  not 
been  detected  (or  overlooked)  the  authors  generally  do  not  pro¬ 
vide  upper  limits.  Yet,  even  so,  the  biases  affecting  our  sources 
selection  do  not  prevent  us  from  analysing  the  outflow  proper¬ 
ties  within  this  sample,  which  is  nevertheless  an  interesting  and 
quite  heterogeneous  collection  of  local  galaxies,  spanning  a  wide 
range  of  physical  properties,  especially  in  terms  of  star  formation 
rate  and  AGN  luminosity. 

All  of  the  plots  shown  in  this  section  are  obtained  by  assum¬ 
ing  a  CO-to-H2  conversion  factor  trco(i-o)=0.8  for  the  molecular 
gas  in  the  outflow  (see  also  Table  [4]  and  the  related  discussion 
in  Section  3.2,  for  the  PdBI  sample,  and  Section  4.2,  for  the 
literature  sample);  nevertheless,  the  use  of  the  alternative  con¬ 
version  factors  (whose  resulting  outflow  rates  and  energetics  are 
also  reported  in  Table]?])  does  not  affect  the  observed  trends  sig¬ 
nificantly.  The  errors  on  the  outflow  radius  (i.e.  ±0.1  dex,  on 
average),  outflow  rate  and  depletion  time- scale  (±0.3  dex),  out¬ 
flow  kinetic  power  (±  0.5  dex)  and  momentum  rate  (±  0.45  dex), 
are  calculated  for  our  4  PdBI  sources  with  detected  molecular 
outflows,  by  using  the  uncertainties  on  Sco  and  Rof  obtained 
from  the  uv  data  (reported  in  Table[2]),  and  by  assuming  a  conser¬ 
vative  average  outflow  velocity  error  of  ±50%,  which  accounts 
for  uncertainties  associated  with  projection  effects.  We  assume 
these  error  estimates  to  be  representative  also  of  the  sample  ex¬ 
tracted  from  the  literature,  for  which  obtaining  reliable  error- 
bars  is  quite  difficult,  given  the  incompleteness  of  the  informa¬ 
tion  provided  in  the  literature.  These  error-bars  do  not  take  into 
account  the  uncertainty  on  the  CO-to-H2  conversion  factor;  how¬ 
ever,  we  note  that  including  it  would  simply  increase  the  upper 
error-bars,  since  we  are  using,  for  all  the  sources,  the  most  con¬ 
servative  value  of  (Tcoq-o)  (in  the  molecular  outflow). 


6.1.  The  Relation  Between  Outflow  Rate,  SFR  and  AGN 
Luminosity. 

The  top  panel  of  Figure  [8]  shows  the  outflow  rate  Mh2,of  as  a 
function  of  the  star  formation  rate,  by  colour-coding  the  different 
symbols  according  to  the  fraction  of  the  total  galaxy  bolometric 
luminosity  that  is  ascribed  to  the  AGN  (i.e.  Lagn/LboO-  Solid 
and  open  circles  indicate  galaxies  optically  classified  as  type 
1  (unobscured)  and  type  2  (obscured)  AGNs,  respectively;  tri¬ 
angles  represent  LINER-type  galaxies  and  stars  indicate  galax¬ 
ies  classified  as  “pure”  starbursts  (HII).  In  this  diagram,  the 
pure  starbursts  and  the  starburst-dominated  galaxies  lie  close 
to  the  relation  SFR  =  Mh2,of>  he.  outflow  mass  loading  factor 
rj  =  Mof/SFR  ~  1.  Conversely,  objects  with  higher  AGN  frac¬ 
tion  depart  from  this  sequence,  and  the  deviation  appears  to  in¬ 
crease  with  LAgn/Lbo1-  This  result  would  suggest  that  the  pres¬ 
ence  of  an  AGN  can  boost  the  outflow  rate  by  a  large  factor, 
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Fig.  8:  Top:  Outflow  mass-loss  rate  as  a  function  of  the  star  for¬ 
mation  rate  for  the  extended  sample  of  galaxies,  whose  general 
properties  and  outflow  characteristics  are  summarised  in  Tables 
[3]  and  [4]  Filled  and  open  circles  represent  respectively  unob¬ 
scured  and  obscured  AGNs,  LINERs  are  plotted  as  upward  tri¬ 
angles  and  “pure”  starburst  galaxies  as  stars  (see  legend  at  the 
top-left  corner  of  the  plot).  Symbols  are  colour-coded  according 
to  the  fraction  of  bolometric  luminosity  attributed  to  the  AGN 
(Lagn/Lboi)-  Outflow  rates  are  derived  by  assuming,  in  the  out¬ 
flow,  the  conservative  CO-to-F^  conversion  factor  trco(i-o)=0.8 
for  all  of  the  sources.  The  black  dashed  line  represents  the  1:1 
correlation  between  SFR  and  outflow  mass-loss  rate.  Bottom: 
This  plot  indicates  a  positive  correlation  between  the  outflow 
mass  loading  factor  (MH2,of/SFR)  and  LAgn/Lb0i  that  emerges 
from  the  diagram  in  the  top  panel.  Notes:  At  the  bottom-right 
of  each  panel  we  show  an  error-bar  representative  for  the  whole 
sample,  which  does  not  take  into  account  the  uncertainty  on  the 
CO-to-H2  conversion  factor  (see  text  in  Section  6  for  a  more  de¬ 
tailed  explanation),  and  corresponds  to  an  average  error  of  ±0.3 
dex. 


which  approaches  two  orders  of  magnitude  in  the  most  powerful 
QSOs  of  our  sample. 

However,  we  note  that  also  some  of  the  SB -dominated 
sources  have  outflow  rates  which  depart  slightly  from  a  1 : 1  cor¬ 
relation  by  having  outflow  mass  loading  factors  77  ~  1-4,  indi¬ 
cating  that  they  are  also  in  a  star  formation-quenching  regime. 
Although  models  of  starburst-driven  outflows  can  predict  mass 
loading  factors  of  up  to  ~2-3,  such  high  values  of  the  outflow 
mass  loading  factor  in  SB -dominated  sources  may  indicate  that 


Article  number,  page  15  of  26 


l°g(LAGN)  (erg  s-1) 
L’AGn/^BqI 


0.00  0.25  0.50  0.75  1.00 

Fig.  9:  Outflow  mass-loss  rate  as  a  function  of  the  AGN  bolo- 
metric  luminosity,  for  the  extended  sample  of  galaxies  (Tables 
|3]and|§.  Symbols  and  colour-coding  are  as  in  Figure  [8]  Error 
bars  as  in  Fig.  [8]  For  the  “pure”  starburst  galaxies  we  set  an 
upper  limit  to  the  AGN  contribution  by  hypothesising  that  all  of 
their  hard  X-ray  (2-10  keV)  unobscured  luminosity  is  produced 
by  an  AGN  (further  explanation  in  Section  5).  The  red  dashed 
relationship  results  from  a  linear  fit  to  the  AGN  host-galaxies,  in 
which  the  upper  limits  have  been  excluded. 


even  in  these  case  an  AGN  can  contribute  in  driving  the  outflow, 
even  if  the  AGN  accounts  only  for  a  minor  fraction  of  the  to¬ 
tal  bolometric  luminosity  of  the  source.  Alternatively,  an  AGN 
may  have  been  present,  but  switched  off  recently.  Indeed,  re¬ 
cent  models  show  that  AGN-driven  massive  molecular  outflows 
can  persist  for  ~  108  yrs  after  the  central  nucleus  has  turned  off 
([King  et  al.|2011  Zubovas  &  King|2012). 

The  tentative  positive  correlation  between  the  mass  loading 
factor  Mh2,of/SFR  and  the  AGN  fraction  of  the  bolometric  lumi¬ 
nosity  of  the  galaxies  is  better  illustrated  by  the  the  bottom  panel 
of  Figure  [8] 

The  relevance  of  the  AGN  in  driving  outflows  is  also  high¬ 
lighted  in  the  plot  of  the  outflow  rate  as  a  function  of  AGN  lumi¬ 
nosity  reported  in  Fig.  [9]  For  the  galaxies  of  our  sample  which 
host  an  AGN,  the  AGN  luminosity  correlates  with  the  outflow 
rate.  Interestingly,  “pure”  starburst  galaxies  are  all  outliers  in 
this  correlation,  possibly  indicating  that  the  feedback  mecha¬ 
nism  in  action  in  these  objects  is  substantially  different  from  the 
AGN-host  galaxies,  even  from  those  AGNs  whose  bolometric 
luminosity  is  dominated  by  the  starburst.  We  fit  the  correlation 
for  the  AGN-host  galaxies,  by  excluding  the  upper  limits,  and 
we  obtain: 


l0g(MH2,0F)  =  (-29.8  ±  3.7)  +  (0.720  ±  0.083)  log(LAGN),  (2) 

resulting  in  the  red  dashed  line  in  Figure  [9]  This  is  the  first 
time  that  a  correlation  between  the  molecular  outflow  rate  and 
the  AGN  luminosity  is  found  from  CO  data,  and  we  interpret 
it  as  a  direct  evidence  for  these  massive  molecular  outflows  be¬ 
ing  mostly  powered  by  the  AGN  and,  therefore,  for  AGN  neg¬ 
ative  feedback  in  action  in  the  AGN-host  galaxies  included  in 
our  study.  However,  I  Zwl  and  IC  5063  are  remarkable  counter 
examples,  deviating  from  the  relation.  Given  the  biases  and  in¬ 
completeness  of  our  sample,  these  two  objects  may  be  the  tip  of 
the  iceberg  of  a  population  of  AGNs  with  low  outflow  rates.  In 


this  scenario,  the  correlation  in  Fig.  [9]  would  be  just  the  upper 
envelope  of  the  entire  distribution  of  AGN  host-galaxies  in  the 
outflow  rate  vs  AGN  luminosity  diagram.  This  possibility  should 
be  investigated  with  extensive  CO  surveys,  targeting  a  complete 
sample  of  galaxies.  Alternatively,  I  Zwl  and  IC  5063  may  rep¬ 
resent  a  transient  phase  in  which  the  molecular  outflow  has  been 
swept  away  or  the  molecules  photo-dissociated  and  ionised,  and 
additional  molecular  clouds  from  the  disk  have  still  to  be  accel¬ 
erated. 

6.2.  Gas  Depletion  Time-Scales 


Fig.  10:  Molecular  gas  depletion  time-scale  due  to  the  molecu¬ 
lar  outflow  (TdeP,oF  =  Mh2,tot/Mh2,of)  versus  AGN  bolometric 
luminosity.  Symbols  and  colour-coding  as  in  Fig.  [8]  The  uncer¬ 
tainties  (±0.3  dex)  on  the  depletion  time- scales  (a  representative 
error-bar  is  shown  at  the  bottom-left  of  the  plot)  are  dominated 
by  the  errors  on  the  outflow  rates. 
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Fig.  11:  Comparison  between  the  depletion  time-scale  due  to 
gas  removal  through  molecular  outflows  with  the  depletion  time- 
scale  due  to  gas  consumption  by  star  formation.  Symbols  and 
colour-coding  as  in  Fig.  [8]  Error-bars  as  in  Fig.  10 


Fig.  [TO]  shows  the  molecular  gas  depletion  time- scale  (due 
to  outflows)  as  a  function  of  the  AGN  luminosity.  The  gas  de- 
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pletion  time-scale  associated  with  the  outflow  is  defined  as  the 
time  required  by  the  outflow  (assuming  that  it  continues  at  the 
same  rate)  to  completely  remove  the  whole  gas  content  in  the 


host  galaxy,  i.e.  rdep5oF  =  MH2,tot/MH2,of.  In  Fig.  10  the  gas 


depletion  time- scales  are  significantly  reduced  by  the  presence 
of  a  luminous  AGN.  We  observe  indeed  an  anti-correlation  be¬ 
tween  these  two  quantities,  as  already  observed  by  Sturm  et  al. 
( 201 1|),  implying  that  the  objects  hosting  more  powerful  AGNs 
are  depleted  of  their  molecular  gas  content  on  shorter  time- 
scales.  Starburst  galaxies  have  outflow  depletion  time- scales  up 
to  several  hundred  million  years,  while  in  galaxies  hosting  pow¬ 
erful  quasars  the  depletion  time- scales  can  be  as  short  as  a  few 
million  years.  We  note  however  that  these  are  time-scales,  which 
would  correspond  to  the  real  depletion  times  only  if  the  outflow 
continues  at  this  rate.  Moreover,  as  already  stressed  at  the  begin¬ 
ning  of  this  section,  our  sample  is  strongly  biased  toward  galax¬ 
ies  and  quasars  which  do  host  massive  molecular  outflows,  and 
therefore  these  results  may  not  hold  for  the  general  population 
of  local  galaxies. 

Fig.  [IT] compares  the  depletion  time- scale  due  to  the  molec¬ 
ular  outflow  with  the  depletion  time- scale  due  to  gas  consump¬ 
tion  by  the  star  formation.  The  line  indicates  the  locus  where  the 
two  time-scales  are  identical.  In  starburst-dominated  galaxies 
the  two  time- scales  are  similar,  while  in  AGN-dominated  galax¬ 
ies  the  outflow  depletion  time- scale  is  much  shorter  than  the  de¬ 
pletion  time-scale  due  to  star  formation. 

It  is  interesting  to  note  that  outflows  can  expel  the  molecular 
gas  on  time-scales  shorter  than  100  Myr,  which  is  the  minimum 
quenching  time- scale  required  to  explain  the  ^-enhancement  in 
massive  spheroidal  galaxies.  Therefore,  if  the  same  mechanism 
is  in  place  at  high  redshift,  it  can  actually  help  to  explain  the 
enhancement  of  tr-elements,  relative  to  iron,  observed  in  local 
elliptical  galaxies.  It  is  also  worth  noting  that  the  outflow  deple¬ 
tion  time-scale,  especially  in  the  powerful  AGNs  of  our  sample, 
can  be  even  shorter  than  10-20  Myr,  which  is  the  quenching 
time  inferred  by  recent  studies  of  post-starburst  galaxies  at  in¬ 
termediate  and  high  redshift  (Christy  Tremonti,  Priv.  Comm.). 
Conversely,  the  depletion  time- scales  associated  with  the  con¬ 
sumption  of  gas  by  star  formation  are  too  long  and  fail  to  meet 
these  conditions  in  most  objects. 


6.3.  Kinetic  Power  of  the  Outflows 

Fig.  [12]  shows  the  kinetic  power  of  the  molecular  outflow  as  a 
function  of  AGN  luminosity.  Theoretical  models  of  AGN  feed¬ 
back  and  cosmological  simulations  predict  a  coupling  efficiency 
between  AGN-driven  outflows  and  AGN  power  of  about  ~5%, 
for  AGN  accreting  close  to  the  Eddington  limit  (which  is  likely 
the  case  for,  at  least,  the  most  luminous  AGNs  in  our  sample). 
This  is  also  the  Pkin/LAGN  fraction  needed  to  explain  the  Mbh-0" 
relation  in  local  galaxies  (e.g.  |King|2010l|Zubovas  &  King|20T2| 


|Lapi  et  al71|2005|>.  Our  observations  of  massive  molecular  out 


flows  in  AGN-host  galaxies,  overall,  appear  to  confirm  this  pre¬ 
diction.  The  solid  grey  line  in  Fig.  [12]  indicates  the  locus  of 
points  having  an  outflow  kinetic  power  that  is  5%  of  the  AGN 
luminosity,  and  galaxies  hosting  powerful  quasars  are  indeed  lo¬ 
cated  close  to  this  value. 

The  best-fit  to  our  data  points  by  excluding  the  upper  limits 
is: 


log(Pkin,oF)  =  (-9.6  ±  6.1)  +  (1.18  ±  0.14)  log(LAGN)  (3) 

and  it  is  indicated  by  the  red  dashed  line  in  Figure  [12]  It  is  in¬ 
teresting  that  low  luminosity  AGNs  seem  to  show  an  efficiency 
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Fig.  12:  Correlation  between  the  kinetic  power  of  the  outflow 
and  the  AGN  bolometric  luminosity.  Symbols  and  colour-coding 
as  in  Fig.  [8]  The  grey  line  represents  the  theoretical  expectation 
of  models  of  AGN  feedback,  for  which  Pk,of  =  5%LAgn-  The 
red  dashed  line  represents  the  linear  fit  to  our  data,  excluding  the 
upper  limits.  The  error-bar  shown  at  the  bottom-right  of  the  plot 
corresponds  to  an  average  error  of  ±0.5  dex. 


41.5  42.0  42.5  43.0  43.5  44.0 
log(PK)sF  (erg  s_1) 

LaGn/FboI 

0.00  0.25  0.50  0.75  1.00 

Fig.  13:  Kinetic  power  of  the  outflow  plotted  as  a  function  of  the 
kinetic  power  of  a  supernova-driven  wind.  Symbols  and  colour¬ 
coding  as  in  Fig.  [8]  Error  bars  as  in  Fig.  [12]  The  black  and 
grey  lines  mark  the  relations  Pk,of  =  Pk,sf,  Pk,of  =  10%  Pk,sf 
andPK,oF  =  1%Pk,sf- 


lower  than  5%.  Likely,  these  AGNs  (especially  the  LINERs,  in¬ 
dicated  with  triangles)  are  accreting  at  a  rate  lower  than  Edding¬ 
ton.  One  should  also  note  that  in  some  of  these  low  luminosity 
AGNs  (e.g.  IC  5063,  see|Morganti  et  al.|  d2013|),  and  possibly 
NGC  1266  and  NGC  6764,  as  suggested  by|Alatalo  et  alJpOTTT) 
and|Leon  et  al.|(|2007]),  respectively)  a  radio-jet  is  thought  to  con¬ 
tribute  to  the  acceleration  of  the  molecular  gas.  Additional  de¬ 
tailed  observations  are  required  to  better  understand  the  outflow 
driving  mechanism  in  these  objects. 

Starburst  galaxy  upper  limits  are  located  above  the  5%  line, 
indicating  that  in  these  objects  a  different  source  of  energy  is  re- 
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Fig.  14:  Outflow  momentum  rate  ((v  Mh2,of))  versus  photon 
momentum  output  of  the  AGN  (LAgn/c).  Symbols  and  colour¬ 
coding  as  in  Fig.  [8]  The  grey  line  shows  the  prediction  of  mod¬ 
els  of  AGN  feedback,  i.e.  Mh2,of  v  ~  20  Lagn/c.  The  error-bar 
shown  at  the  bottom-right  of  the  plot  corresponds  to  an  average 
error  of  ±0.45  dex. 


quired,  most  likely  provided  by  SN  ejecta  and  radiation  pressure 
from  the  young  stars. 

The  kinetic  power  of  the  outflow  is  compared  in  Figure  [13] 
with  the  kinetic  power  injected  by  supernovae,  as  inferred  from 
the  SFR,  following  |Veilleux  et  al.|  (|2005|)  (see  Table  [3]  and  rel¬ 
evant  explanation  in  Section  5).  Fig.  [13]  shows  that  the  out¬ 
flow  kinetic  power  achieved  in  the  “pure”  starburst  galaxies  and 
in  some  of  the  starburst-dominated  objects  is  compatible  with  a 
supernova-driven  wind,  with  a  coupling  efficiency  of  a  few  per¬ 
cent  up  to  a  few  10  percent.  Conversely,  for  most  of  the  heavily 
AGN-dominated  sources,  it  is  clear  that  the  additional  contri¬ 
bution  of  the  AGN  is  needed  to  produce  the  observed  outflow 
energetics. 


6.4.  Momentum  Rate  of  the  Outflows 


The  momentum  rate  provides  an  additional  important  indicator 
of  the  nature  of  the  outflow  and  an  important  test  for  models. 
In  models  in  which  the  outflow  is  generated  by  a  nuclear  AGN- 
driven  wind,  the  momentum  rate  transferred  by  the  AGN  pho¬ 
tons  to  the  surrounding  medium  is  given  by  the  average  number 
of  scattering  by  each  photon.  Some  of  these  models  predict, 
for  AGNs  accreting  close  to  the  Eddington  limit,  momentum 
fluxes  of  the  order  of  ~  20  LAgn/c  (e.g.  |Zubovas  &  King|2012[ 
Faucher-Giguere  &  Quataert  2012). 

Very  interestingly,  we  find  that  the  “momentum  boost”,  i.e. 
the  ratio  of  vMh2,of  to  the  AGN  radiative  momentum  output 
Lagn/c,  ranges  from  ~10  to  ~50  in  the  galaxies  where  the 
AGN  contributes  more  than  10%  to  the  total  bolometric  lumi¬ 
nosity  (Table  [4]).  In  particular,  most  of  our  sources  (except 
the  “pure”  starbursts)  do  follow,  within  the  errors,  the  relation 


vMh2,of  ~  20  Lagn/c  (Fig.  14).  This  finding  both  supports 


the  AGN  energy-driven  nature  of  these  outflows  and  the  AGN 
feedback  models  that  have  been  proposed  so  far,  in  which  a 
fast  and  highly  ionised  wind,  arising  from  the  nuclear  regions  of 
the  AGN,  creates  a  shock  wave  that  propagates  into  the  ISM  of 
the  galaxy  ([Zubovas  &  King|2012[  Faucher-Giguere  &  Quataert 
2012  ).  We  further  explore  this  hypothesis  in  Fig.  15  by  repro¬ 


Fig.  15:  Outflow  momentum  “boost"  with  respect  to  the  AGN 
radiative  momentum  output,  i.e.  (v  Mh2, of) /(Lagn/c),  as  a  func¬ 
tion  of  the  outflow  velocity,  following  Fig.  5  of  Faucher-Giguere 
&  Quataert 1 12012|).  Symbols  and  colour-coding  as  in  Fig.  [8]  The 


uncertainties  on  the  momentum  boost  are  dominated  by  the  er¬ 
rors  on  the  outflow  momentum  rates.  The  upper  limits  on  the 
outflow  detection  are  excluded  from  this  plot,  as  well  as  the 
“pure"  starbursts. 


ducing  the  diagram  originally  proposed  by  [Faucher-Giguere  &] 
Quataert  (  2012  )  (Figure  5  in  their  paper),  in  which  the  momen¬ 
tum  “boost”  of  the  outflows  is  plotted  as_a  function  of  the  out¬ 
flow  velocity.  According  to  [Faucher-Giguere  &  Quataert  ( 2012 ), 
under  the  assumption  that  the  observed  molecular  outflows  are 
energy-conserving,  the  momentum  boost  Mh2, of /(Lagn/c)  and 
the  molecular  outflow  velocity  vof  are  linked  to  the  velocity  of 
the  initial  nuclear  and  fast  wind,  Vin,  by  the  relationship: 


Mh2,OF  1  V^ 

(Lagv/c)  2  vqf 


(4) 


This  relationship,  computed  for  different  values  of  Vin,  is  rep¬ 
resented  by  the  grey  dashed  lines  in  Fig.  |T5]  It  is  clear  from 
Fig.  [I5]that  we  do  not  observe  the  correlation  found  by|Faucher^ 
Giguere  &  Quataert  (2012);  however,  the  trend  shown  in  their 
paper  is  mainly  driven  by  the  two  sources  with  outflow  veloc¬ 
ities  of  ~5000  km  s-1,  in  which  the  powerful  winds  have  been 
observed  in  ionised  gas  (FeLoBALs).  Due  to  the  narrow  range  of 
velocities  sampled  by  our  molecular  outflows,  we  probably  can¬ 
not  test  the  correlation  of  |Faucher-Giguere  &  Quataert]  (|20 12]), 
but  we  can  still  comment  on  the  consistency  of  their  model  with 
our  observations.  Fig [15]  illustrates  that  nuclear  ionised  winds 
with  Vin  e  (0.01c,  0.1c)  can  in  principle  generate  the  molecular 
outflows  observed  in  the  AGN-dominated  sources. 

However,  we  should  still  keep  in  mind  the  bias  affecting  our 
sample,  which  has  been  selected  to  mostly  include  objects  with 
known  (powerful  outflows).  Within  this  context  it  is  useful  to 
note  that  the  upper  limit  given  by  I  Zw  1  points  towards  a  much 
lower  “momentum  boost”  in  the  plot  in  Fig.  PM]  If  this  object 


is  representative  of  a  wider  population  (under-represented  in  our 
work  because  of  our  selection  biases),  then  it  may  hint  at  other 
outflow  acceleration  mechanisms  at  work  in  AGNs,  such  as  di¬ 
rect  radiation  pressure  onto  the  dusty  clouds  in  the  host  galaxy 
(e.g.  |Fabian||2012|).  Moreover,  if  the  distribution  of  the  clouds 
in  the  outflow  is  different  with  respect  to  what  assumed  by  us  in 
Sect.  3.2  (i.e.  continuous  versus  explosive  scenarios)  then  the 
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Fig.  16:  Outflow  momentum  rate  versus  total  photon  momentum 
output  of  the  galaxy  (LBoi/c).  Symbols  and  colour-coding  as  in 
Fig.  [8]  Error  bars  are  as  in  Fig.  [14]  The  grey  line  shows  the 
relation  Mh2,of  v  ~  20  Lboi/c. 


log(v0F)  (km  s"1) 

^AGn/^BoI 

0.00  0.25  0.50  0.75  1.00 

Fig.  17:  Outflow  mass  loading  factor  (Mh2,of/SFR)  plotted  as  a 
function  of  the  outflow  velocity.  Only  the  objects  with  an  AGN 
contribution  to  the  total  bolometric  luminosity  of  <50%  are  plot¬ 
ted.  Symbols  and  colour-coding  as  in  Fig.  [8]  Error  bars  as  in  Fig. 
[8]  The  red  dashed  line  represents  the  linear  fit  to  the  five  “pure" 
starbursts  galaxies  of  our  sample  (M82,  NGC  3256,  NGC  3628, 
NGC  253,  NGC  2146). 


|2012|).  In  this  final  section  we  investigate  the  outflow  driving 
mechanism  in  the  starburst  galaxies  included  in  our  sample  by 
exploiting  our  results. 

It  has  already  been  shown  in  Fig. 


13  that  the  kinetic  power 


of  the  outflow  in  starburst  galaxies  corresponds  to  only  about 
1%  of  the  kinetic  power  produced  by  the  SNe  associated  with 
the  starburst.  Therefore,  if  supernovae  power  the  outflows,  the 
coupling  efficiency  between  supernova  shocks  and  ISM  must  be 
quite  low. 

Additional  information  is  given  by  the  investigation  of  the 
outflow  momentum  rates.  In  particular,  it  is  interesting  to  com¬ 
pare  the  outflow  momentum  rates  with  the  total  radiative  mo¬ 
mentum  output  of  these  galaxies,  i.e.  with  LBoi/c.  For  the  heav¬ 
ily  AGN-dominated  sources,  for  which  Lb0i  ~  Lagn,  such  a  re¬ 
lation,  shown  in  Fig.  [16]  naturally  resembles  that  of  Fig.  [14] 
showing  momentum  “boosts"  of  up  to  a  factor  of  20.  However, 


for  the  other  AGN  hosts  (with  lower  AGN  contribution),  Fig.  16 


demonstrates  that,  when  comparing  the  outflow  momentum  rate 
with  Lsoi/c,  i.e.  when  including  the  contribution  of  the  starburst 
luminosity  in  the  photon  momentum  computation,  “momentum 
boosts”  greater  than  ~20  are  no  longer  required  to  explain  our 
observations.  Instead,  a  significant  contribution  to  the  momen¬ 
tum  rates  of  these  molecular  outflows  may  come  from  momen¬ 
tum  deposition  into  the  ISM  from  star  formation.  In  short,  in 
these  objects  probably  the  radiation  pressure  from  the  AGN  and 
from  starbursts  synergically  combine  to  drive  the  outflow. 

For  the  five  “pure”  starburst  galaxies,  (along  with  the  low 
luminosity  AGN  NGC  6764,  which  however  has  an  AGN  con¬ 
tribution  of  only  ~2%)  Fig.  [16]  shows  that  the  momentum  rate 


is  comparable  to  LBoi/c,  suggesting  that  (single- scattering)  radi¬ 
ation  pressure  onto  the  dusty  clouds  of  the  ISM  may  contribute 
substantially  to,  or  even  dominate,  the  driving  mechanism  of  out¬ 
flows  in  starbursts. 

The  latter  scenario  is  further  supported  by  Fig.  [IT]  which 
shows  the  outflow  mass  loading  factor  (MH2,of/SFR)  as  a 
function  of  outflow  velocity.  In  energy-driven  outflows  the 
mass  loading  factor  should  be  proportional  to  v^,  while  in 
momentum-driven  outflows  the  mass  loading  factor  should  be 


proportional  to  vQI(e.g.  Murray  et  al 


2005).  The  two  dashed 


gray  lines  in  Fig.  [IT]  show  a  tit  to  the  starburst  galaxies  with  the 
slopes  fixed  to  -1  and  -2.  Clearly  the  fit  with  slope  -2,  which 
would  be  compliant  with  the  energy-driven  scenario,  is  incon¬ 
sistent  with  the  starburst  data.  Instead,  the  fit  with  slope  -1, 
expected  from  momentum-driven  models,  is  in  nice  agreement 
with  the  starburst  data.  The  red  dashed  line  shows  the  best  fit 
to  the  starburst  data  obtained  by  leaving  the  slope  free,  which 
gives: 


momentum  rate  may  be  lower  by  a  factor  of  three,  hence  possi¬ 
bly  favouring  the  latter  models. 


6.5.  The  Driving  Mechanism  in  the  “Pure”  Starbursts 


In  all  of  the  previous  sections,  we  have  seen  that  “pure”  star- 
burst  galaxies  deviate  from  the  outflow  relations  found  for  galax¬ 
ies  hosting  AGNs,  indicating  that  in  these  systems  the  out¬ 
flow  driving  mechanism  is  likely  different.  Two  main  outflow 
mechanisms  have  been  proposed  for  starburst  galaxies:  “energy- 
driven”,  commonly  associated  with  the  energy  injected  into  the 
ISM  by  SNe  explosions,  and  “momentum  driven”,  generally  as¬ 
sociated  with  radiation  pressure  from  the  young  stellar  popula¬ 
tion  on  the  dusty  clouds  ([Murray  et  al.[2005,  2011;  Thompson 


et  al. ||2005l  |Oppenheimer  &  Dave||2006|  |2008[  |Hopkins  et  al. 


log(MH2,0F/SFR)  =  (1.9  ±  1.0)  +  (-1.0  ±  0.5)  log(vOF)  (5) 

As  shown  by  Eq.  (5),  the  best  fit  to  the  pure  starbursts  is 
fully  consistent  with  the  slope  -1,  while  it  is  only  marginally 
consistent  (2 cr)  with  slope  -2.  In  summary,  our  results,  although 
based  on  a  very  modest  sample  size,  suggest  that  in  starburst 
galaxies  outflows  are  mostly  momentum-driven,  at  least  for  the 
molecular  component  (which  generally  dominates  the  outflow 
mass  budget). 

7.  Summary  and  Conclusions 

We  have  investigated  the  properties  of  massive  molecular  out¬ 
flows  in  a  sample  of  19  nearby  galaxies,  spanning  a  wide  range 
of  AGN  and  starburst  activities.  For  seven  of  these  sources,  we 
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have  presented  new  PdBI  observations  of  the  CO(l-O)  emission 
line;  in  four  of  them,  we  have  detected,  for  the  first  time,  CO 
broad  wings  tracing  molecular  outflows  extended  on  kpc  scales, 
which  indicates  that  these  outflows  are  affecting  their  host  galax¬ 
ies  on  large  scales.  In  other  two  sources  we  have  found  some 
marginal  evidence  for  the  presence  of  broad  CO  wings,  but  ad¬ 
ditional  observations  are  required  to  confirm  the  presence  and 
the  nature  of  this  high  velocity  gas.  We  have  measured  very 
large  outflow  mass-loss  rates  of  several  times  100  M©  yr_1,  sim¬ 
ilarly  to^heprototypical^iantniolecularwind  discovered  in  Mrk 
231  ([Feruglio  et  al.|2010t|Fischer  et  al.|2010|).  With  the  aim  to 
shed  light  on  the  possible  correlations  between  outflow  charac¬ 
teristics  and  intrinsic  properties  of  a  larger  sample  of  objects,  we 
have  included  in  our  study  12  additional  local  galaxies,  in  which 
outflows  of  molecular  gas  have  been  observed.  For  homogene¬ 
ity,  we  have  restricted  our  investigation  to  molecular  outflows  re¬ 
vealed  by  observations  of  high  velocity  components  of  the  CO(l- 
0)  or  CO(2-l)  emission  lines.  We  emphasize  that  the  sample 
employed  for  our  analysis  (both  our  new  sample  and  the  sam¬ 
ple  taken  from  the  literature)  is  clearly  biased  towards  galaxies 
for  which  a  molecular  outflow  has  been  detected.  By  comparing 
the  intrinsic  characteristics  of  this  sample  of  nearby  galaxies  and 
AGN  hosts  with  their  outflow  properties  and  energetics,  we  have 
found  that: 


-  Starbursts  are  effective  at  powering  massive  molecular  out¬ 
flows,  but  the  presence  of  a  powerful  AGN,  when  it 
dominates  the  energetic  output  of  the  galaxy,  can  signifi¬ 
cantly  boost  the  outflow  rates.  More  specifically,  while  in 
starburst-dominated  galaxies  the  outflow  mass  loading  fac¬ 
tor,  ( rf  =  Mh2,of/SFR),  is  close  to  unity,  with  maximum  val¬ 
ues  of  2-4,  the  mass  loading  factor  rj  is  found  to  increase 
with  the  fraction  of  total  bolometric  luminosity  that  is  as¬ 
sociated  with  the  AGN,  approaching  values  of  -100  in  the 
most  powerful  QSOs. 

-  We  have  found  a  correlation  between  outflow  rate  and  AGN 
luminosity  for  the  AGN  host  galaxies  of  our  extended  sam¬ 
ple.  However,  this  correlation  may  be  tracing  the  upper  enve¬ 
lope  of  a  broader  distribution.  Indeed  our  sample  is  certainly 
biased  towards  objects  whose  presence  of  a  molecular  out¬ 
flow  was  known  a  priori.  A  completely  unbiased  survey  may 
deliver  more  sources  with  lower  outflow  rates,  placed  below 
the  relation  found  by  us.  The  presence  of  two  objects,  in  our 
sample,  hosting  powerful  AGNs  and  placed  below  the  rela¬ 
tion,  hints  at  this  scenario. 

-  We  have  found  that  these  AGN-driven  massive  molecular 
outflows  are  capable  of  exhausting  the  cold  gas  reservoir  of 
the  host  galaxy  and,  hence,  of  quenching  star  formation,  on 
time-scales  as  short  as  a  few  times  106  yrs,  for  the  most  lumi¬ 
nous  AGNs.  In  particular,  we  have  found  an  anti-correlation 
between  depletion  time- scale  and  AGN  power,  confirming 
previous  studies  based  on  the  far-IR  OH  absorption  lines 
(Sturm  et  al.|2011|).  In  AGN-hosts,  the  depletion  time-scale 
due  to  the  molecular  outflow  is  much  shorter  than  the  deple¬ 
tion  time-scale  of  gas  consumption  due  to  star  formation. 

-  A  comparison  between  kinetic  power  of  the  outflow  and 
AGN  power  has  shown  that  these  two  quantities  are  corre¬ 
lated,  and  that  their  ratio,  especially  for  the  most  powerful 
AGNs  (likely  accreting  close  to  the  Eddington  limit),  is  very 
close  to  the  value  of  5%  predicted  by  theoretical  models, 
commonly  adopted  in  simulations  of  AGN  negative  feed¬ 
back,  and  needed  to  match  the  MBh  -  cr  relationship.  Fur¬ 
ther  confirmation  to  the  hypothesis  of  AGN-driven  outflows 
is  provided  by  the  finding  that  the  energetics  of  the  molec¬ 


ular  outflows  is  totally  inconsistent  with  supernovae-driven 
winds  in  AGN-dominated  sources. 

-  We  find  that  outflow  momentum  rates  of  ~  20LAgn/c  are 
common  among  the  AGN-host  galaxies  of  our  extended  sam¬ 
ple.  This  likely  indicates  that  the  molecular  outflows  that 
we  observe  in  these  sources  are  energy-conserving,  and  they 
may  have  been  generated  by  fast  nuclear  winds  that  subse¬ 
quently  impact  the  interstellar  medium  of  the  host  galaxy. 
However,  other  driving  mechanisms  (such  as  direct  acceler¬ 
ation  of  dusty  clouds  through  radiation  pressure),  may  also 
contribute,  especially  for  outflows  with  lower  “momentum 
boosts”,  whose  fraction  in  our  sample  may  be  heavily  under¬ 
represented  due  to  selection  biases.  Moreover  one  should 
take  into  account  that  the  distribution  of  clouds  in  the  out¬ 
flow  may  be  different  from  what  assumed  by  us,  and  this 
may  reduce  the  inferred  momentum  rate  by  a  factor  of  about 
three. 

-  Interestingly,  in  the  starburst-dominated  galaxies  (i.e.  those 
where  an  AGN  is  present,  but  it  only  contributes  less  than 
10%  to  the  total  bolometric  luminosity  of  the  galaxy),  the 
power  source  of  the  molecular  outflow  is  ambiguous.  On  the 
one  hand,  the  outflow  rates  are  not  much  larger  than  their 
SFRs.  On  the  other  hand,  in  terms  of  outflow  energetics  and 
momentum  rates,  many  of  them  follow  the  same  relation  as 
AGN-dominated  galaxies,  strongly  suggesting  that  even  in 
these  SB -dominated  objects  the  AGN  can  be  the  main  driv¬ 
ing  source  of  the  outflow.  This  suggests  that  even  if  the 
AGN  is  not  the  main  source  of  energy,  it  can  be  more  ef¬ 
fective  than  the  starburst  in  driving  outflows.  Alternatively,  a 
more  powerful  AGN  may  have  been  present  in  the  past  and 
it  may  have  been  responsible  for  driving  most  of  the  outflow, 
which  can  continue  for  a  few  hundred  million  years  since  the 
AGN  has  faded,  as  predicted  by  some  models.  This  would 
also  explain  the  observations  of  high  velocity  (i.e.  v~1000 
km  s-1)  molecular  outflows  in  stacked  spectra  of  galaxies 
which  do  not  show  a  clear  AGN  contribution  (  Chung  et  al. 
201 1|),  as  well  as  the  detection  of  very  fast  atomic  outflows 
(v  <  1000  km  s-1)  in  galaxies  at  z~0.6  showing  no  evidence 
of  strong  AGN  activity  ( Diamond- Stanic  et  al.  |2012|):  in 
these  cases  the  AGN  may  not  dominate  the  total  luminosity 
(and  its  detection  may  be  made  difficult  by  the  dilution  from 
the  starburst  in  the  host  galaxy),  but  it  may  still  dominate 
the  outflow  mechanism.  AGN  variability  can  also  account 
for  the  lack  of  AGN  detection  in  these  sources.  Moreover,  a 
bias  against  detecting  powerful  AGNs,  even  based  on  mid-IR 
or  X-ray  diagnostics,  in  outflows  oriented  perpendicular  to 
our  line  of  sight  or  during  the  early  stages  of  outflow  (when 
the  AGN  is  still  mostly  embedded  in  high  column  densities), 
may  affect  the  classification  of  some  objects. 

-  For  “pure”  starburst  galaxies  (i.e.  galaxies  which  do  not 
host  an  AGN)  our  data  suggest  that  outflows  (at  least  their 
molecular  phase)  are  probably  momentum-driven.  This  is 
supported  by  the  finding  that  the  outflow  mass  loading  fac¬ 
tor  is  proportional  to  v_1,  as  expected  by  momentum-driven 
models,  and  by  the  finding  that  their  momentum  rate  is  com¬ 
parable  to  Lboi/c,  suggesting  that  (single  scattering)  radiation 
pressure  onto  dusty  clouds  is  responsible  for  expelling  these 
clouds. 

In  conclusion,  we  have  studied  massive  molecular  outflows  in 
galaxies  spanning  a  wide  range  of  starburst  and  AGN  activ¬ 
ities.  In  most  cases,  these  outflows  are  powerful  enough  to 
quench  the  star  formation  in  the  host  galaxy,  by  removing  the 
cold  molecular  gas  out  of  which  stars  form.  Our  results  suggest 


Article  number,  page  20  of  26 


C.  Cicone  et  al.:  Massive  Molecular  Outflows  and  Evidence  for  AGN  Feedback  from  CO  Observations 


that  such  negative  feedback  on  star  formation  is  most  efficient 
and  dramatic  in  powerful  AGNs,  whose  presence  can  signifi¬ 
cantly  boost  outflow  rates,  kinetic  powers  and  momentum  rates 
of  these  molecular  winds  and  reduce  the  depletion  time- scale 
of  molecular  gas.  Starbursts  are  also  effective  in  driving  mas¬ 
sive  outflows  (likely  through  radiation  pressure),  although  not 
as  much  as  AGNs.  However,  one  shall  take  into  account  that 
star  formation  likely  lasts  for  a  longer  period  than  the  AGN  ac¬ 
tivity.  As  a  consequence,  while  AGN  feedback  has  likely  a  dra¬ 
matic  effect  in  rapidly  clearing  massive  galaxies  of  their  gas  con¬ 
tent  on  short  time- scales,  the  integrated  effect  of  starburst-driven 
winds  may  be  more  important  on  longer  time- scales,  especially 
for  lower  mass  galaxies.  Further  observations,  carried  on  with 
the  new  generation  of  millimiter-wave  interferometers  (ALMA, 
NOEMA)  will  allow  us  to  extend  this  investigation  to  a  much 
larger  sample  of  local  galaxies,  and,  most  importantly,  overcome 
the  biases  that  are  affecting  the  current  sample. 
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Appendix  A:  Additional  Information 

We  report  in  this  appendix  additional  information  about  the 
seven  galaxies  that  we  observed  with  the  PdBI  in  their  CO(l-O) 
transition.  We  further  discuss  our  PdBI  CO(l-O)  observations, 
by  showing  the  results  of  the  fits  to  the  uv  plots  of  the  CO  wings, 
which  provide  an  estimate  of  the  size  of  the  molecular  outflow 
(when  detected).  We  also  present  the  integrated  map,  the  corre¬ 
sponding  uv  plot  and  the  velocity  field  of  the  narrow  core  CO(l- 
0)  emission,  which  traces  the  bulk  of  the  molecular  gas  in  rota¬ 
tion  in  the  galaxy  disk/ring. 

Appendix  A.  1 :  IRAS  F08572+3915 


Fig.  A.l:  Real  part  of  visibilities  as  a  function  of  the  uv  dis¬ 
tance  for  the  blue  (left)  and  red  (right)  CO(l-O)  wings  of  IRAS 
F08572+3915,  binned  in  baseline  steps  of  20m.  The  red  curve 
is  the  best  fit  with  a  circular  Gaussian  model.  Note  that  the  uv 
tables  of  the  two  CO(l-O)  wings  have  been  re-centred  to  their 
emission  centroid  positions  (obtained  from  their  maps  shown  in 
Fig.  [T])  prior  to  plotting  the  uv  visibilities  in  order  to  correctly 
estimate  their  size. 


Fig.  A.2:  Panel  (a):  IRAM  PdBI  continuum- subtracted  map  and 
of  the  CO(l-O)  narrow  core  emission  in  IRAS  F08572+3915. 
Contour  levels  correspond  to  5 <x.  The  inset  (b)  shows  the  first 
moment  map  obtained  within  ve(-150,  150)  km  s-1.  Contours 
correspond  to  3  km  s-1.  The  size  of  the  inset  is  6x7  arcsec. 
Panel  (c):  uv  plot  for  the  narrow  core,  with  visibilities  binned  in 
uv  radius  steps  of  20m. 


The  NW  galaxy  hosts  powerful,  multi-phase  and  large-scale 
outflows,  characterised  by  very  high  gas  velocities  in  each  com¬ 
ponent.  The  ionised  and  neutral  winds  have  been  resolved  at 
-kpc  scales  by  Rupke  &  Veilleux  (  2013b),  based  on  several  op¬ 
tical  emission  lines  and  NalD  absorption,  and  they  appear  to  be 
partially  overlapping.  The  ionised  outflow  shows  blue-shifted 
velocities  of  up  to  -3350  km  s-1,  while  the  maximum  blue- 
shifted  velocity  detected  in  NalD  absorption  is  -1150  km  s-1. 
Herschel-PACS  observations  of  IRAS  F08572+3915  revealed 
the  presence  of  a  powerful  molecular  outflow,  with  OH  79  and 


\\9jjm  blue- shifted  velocities  of  up  to  -1100  km  s  1  ([Sturm 


|et  al.||20lT]  | Veilleux  et  al.||2013|.  Interestingly,  recent  Keck 
OSIRIS  observations  of  the  near-IR  H2  roto- vibrational  transi¬ 
tions  show  that  warm  molecular  gas  is  also  outflowing  with  ve¬ 
locities  strikingly  similar  to  those  observed  in  the  cold  molec¬ 
ular  phase  ([Rupke  &  Veilleux]|2013a|).  This  would  suggest  that 
molecular  clouds  may  be  affected  by  shocks  during  the  acceler¬ 
ation  process  in  this  source. 

In  IRAS  F08572+3915,  both  CO(l-O)  wings  are  detected  at 
high  significance  (>  10<x,  see  Fig.  [TJ,  which  allows  us  to  inves¬ 
tigate  their  spatial  extent  separately.  For  this  purpose,  we  fit  the 
visibility  as  a  function  of  the  uv  radius  (the  plot  of  the  real  part  of 
the  visibility  vs  uv  radius  is  shown  in  Figure  [AT])  with  a  circular 
Gaussian  model  and  we  obtain  a  size  (FWHM)  of  1.36  ±  0.45 
kpc  for  the  blue  wing,  and  1.91  ±  0.50  kpc  for  the  red  wing.  The 
value  reported  in  Table [2] is  the  average  of  the  two. 

In  Fig.  |A.2|we  present  the  map  of  the  CO(l-O)  narrow  core 
integrated  within  ±  250  km  s_1,  the  corresponding  uv  plot  and 
the  velocity  field  of  this  central  concentration  of  molecular  gas. 
The  CO(l-O)  narrow  emission  shows  a  non  compact  structure, 
which  is  also  clearly  appreciable  in  the  uv  visibility  data.  Multi¬ 
ple  blobs  of  molecular  gas  at  the  systemic  velocity  are  detected, 
one  of  which  may  be  part  of  the  tidal  bridge  between  the  NW  and 
the  SE  galaxy  of  the  merger.  The  SE  companion  is  discernible 
in  optical  and  near-IR  images,  but  it  has  never  been  detected  in 
molecular  gas  (Evans  et  al.  2002}  iPapadopoulos  et  al.|12012a). 
The  first  moment  map  shows  that  the  molecular  gas  follows  a 
regular  rotation  pattern  oriented  east- west,  with  an  approximate 
rotation  major  axis  at  a  PA=-80  deg.  The  uv  plot  cannot  be  mod¬ 
elled  with  a  single  Gaussian  but  probably  requires  a  more  com¬ 
plex  function. 


Appendix  A.2:  IRAS  FI  0565+2448 


The  Sy2-ULIRG  IRAS  F08572+3915  is  an  interacting  sys¬ 
tem  composed  by  two  nuclei  (NW  and  SE)  separated  by  a  pro¬ 
jected  distance  of  5.4  arcsec  (6.1  kpc);  only  the  north-west  (NW) 
component,  which  is  the  most  luminous  one  and  hosts  a  buried 
QSO,  has  been  detected  in  CO(l-O)  interferometric  maps  ([Evans | 
|et  al.||2002|).  [Papadopoulos  et  ak](|2012b|)  measure  a  very  large 
L^0(6-5)/Lco(i-o)  rati°  the  NW  nucleus,  indicating  very  high 
ISM  excitation,  among  the  most  extreme  in  their  ULIRG  sample. 


Fig.  A.3:  Real  part  of  visibilities  as  a  function  of  the  uv  dis¬ 
tance  for  the  blue  (left)  and  red  (right)  CO(l-O)  wings  of  IRAS 
FI 05 65 +2448,  binned  in  baseline  steps  of  20m.  The  red  curve 
is  the  best  fit  with  a  circular  Gaussian  model.  Note  that  the  uv 
tables  of  the  two  CO(l-O)  wings  have  been  re-centered  to  their 
emission  centroid  positions  (obtained  from  their  maps  shown  in 
Fig.  0  prior  to  plotting  the  uv  visibilities  in  order  to  correctly 
estimate  their  size. 
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Fig.  A.4:  Panel  (a):  IRAM  PdBI  continuum- subtracted  map  of 
the  CO(l-O)  narrow  core  emission  in  IRAS  F10565+2448.  Con¬ 
tours  correspond  to  50<x.  The  inset  (b)  shows  the  first  moment 
map  obtained  within  ve(-200,  200)  km  s-1.  Contours  correspond 
to  3  km  s-1.  The  size  of  this  inset  is  10  x  13  arcsec.  Panel  (c): 
uv  plot  for  the  narrow  core,  with  visibilities  binned  in  uv  radius 
steps  of  5m,  fitted  with  a  power-law  function  model  of  the  form 
oc  r-3  (red  points). 


The  second  Sy2-ULIRG  of  our  sample,  IRAS  F10565+2448, 
is  a  possible  triple  merger  galaxy  system,  whose  emission  is 
largely  dominated  by  the  westernmost  member.  This  is  also 
the  only  merger  component  that  has  been  detected  in  CO(l-O) 
emission  by  Downes  &  Solomon]  (|1998|)  with  the  PdBI  and, 
later  on,  by  |Chung  et  al.  ~(|2009|)  with  the  FCRAO  14m  tele¬ 
scope.  As  for  IRAS  F08572+3915,  the  optical  spectrum  of  IRAS 
F10565+2448  shows  evidence  of  kpc-scale  ionised  and  neutral 
winds,  although  with  outflowing  velocities  much  more  modest 
than  IRAS  F08572+3915  (|Rupke  &  Veilleux|2013b|. 


The  circular  Gaussian  fit  to  the  uv  plots  reported  in  Fig.  A.3 


binned  in  baseline  steps  of  20m,  provides  sizes  (FWHM)  of  2.15 
±  0.32  kpc  and  2.22  ±  0.30  kpc  for  the  blue  and  red-shifted 
CO(l-O)  wings,  respectively. 

The  emission  from  the  CO(l-O)  narrow  core  (i.e.  within  ± 
250  km  s-1),  reported  in  Fig.  A.4i,  appears  more  compact  than 
the  wings,  although  we  find  evidence,  in  the  channel  maps,  for 
a  plume  of  CO(l-O)  emission  extending  to  the  south-west  both 
at  the  systemic  velocity  and  at  slightly  blue- shifted  velocities  of 
about  -150  km  s-1.  Such  plume  is  also  clearly  visible  in  the 
inset  (b)  of  Fig.  |A.4|  which  shows  the  velocity  field  within  the 
velocity  range  of  ve(-200,  200)  km  s-1.  The  major  axis  of  the 
rotation  of  the  molecular  gas  is  oriented  with  a  PA^-60  deg.  The 
plot  of  the  visibility  as  a  function  of  the  uv  radius  proves  that  the 
narrow  component  is  very  well  resolved  in  our  data,  and  its  emis¬ 
sion  can  be  fitted  with  a  power  law  of  the  form  ex  r-3,  resulting 
in  a  FWHM  size  of  0.781  ±  0.078  kpc. 


Appendix  A.3:  IRAS  23365+3604 


Tidal  tails  revealed  by  near-IR  imaging  of  the  LINER-type 
ULIRG  IRAS  23365+3604  indicate  that  this  system  is  in  a  later 
merger  stage,  and  no  evidence  for  a  double  nucleus  has  been 
found  so  far  ([Downes  &  Solomon]  1998|).  Interestingly,  based  on 
a  study  of  the  physical  conditions  of  the  dense  molecular  gas,|Pa-| 
Ipadopoulos  et  al.|([2012a|)  suggest  an  AGN  influence  on  the  ISM 
of  this  source.  IRAS  23365+3604  is  one  of  the  most  spectac¬ 
ular  cases  of  OH  P-Cygni  profiles  detected  by  Herschel-PACS, 
with  blue-shifted  velocities  of  up  to  1300  km  s-1  ( Veilleux  et  al. 

mn). 


In  the  case  of  IRAS  23365+3604,  to  improve  the  signal-to- 
noise  in  order  to  estimate  the  size  of  the  outflow,  we  combine 
the  two  wings  together:  the  corresponding  map  and  uv  plot  are 


100  150 


Fig.  A.5:  Left:  IRAM  PdBI  cleaned  map  of  the  CO(l-O)  blue 
and  red- shifted  wings  of  IRAS  23365+3604,  integrated  within 
(+  300,  +  600)  km  s-1  and  combined  together.  Contour  lev¬ 
els  correspond  to  ler  (  ler  rms  level  is  0.15  mJy  beam-1).  Right: 
Real  part  of  visibilities  as  a  function  of  the  uv  radius  for  the  com¬ 
bined  CO(l-O)  wings;  visibilities  are  binned  in  baseline  steps  of 
50m.  The  red  curve  is  the  best  fit  with  a  circular  Gaussian  model. 


Fig.  A.6:  Panel  (a):  IRAM  PdBI  continuum- subtracted  map  of 
the  CO(l-O)  narrow  core  emission  in  IRAS  23365+3604.  Con¬ 
tours  correspond  to  30<x  (ler  rms  level  is  0.15  mJy  beam-1).  The 
inset  (b)  shows  the  first  moment  map  obtained  within  ve(-200, 
200)  km  s-1.  The  size  of  the  inset  is  8  x  8  arcsec.  Panel  (c): 
uv  plot  for  the  narrow  core,  with  visibilities  binned  in  uv  radius 
steps  of  10m,  fitted  with  a  circular  Gaussian  model. 


shown  in  Figure  |A3|  A  circular  Gaussian  fit  to  the  uv  visibility 
data,  binned  in  baseline  steps  of  50m,  gives  a  FWHM  of  2.45  + 
0.70  kpc,  which  we  use  as  an  estimate  of  the  size  (diameter)  of 
the  molecular  outflow. 

In  Fig.  |A.6fr  we  also  show  the  cleaned  map  of  the  CO(l- 
0)  narrow  core,  integrated  within  +  250  km  s-1;  with  a  flux 
of  Sco,core  “  (45.24  +  0.48)  Jy  km  s-1,  this  component  con¬ 
tributes  more  than  90%  to  the  total  CO(l-O)  emission  of  IRAS 
23365+3604.  Our  observations  resolve  reasonably  well  this  core 
emission,  and  a  circular  Gaussian  fit  to  the  uv  visibilities,  binned 
in  steps  of  10m  of  uv  distance,  results  in  a  FWHM  of  3.174 
+  0.073  kpc.  We  note  in  the  channel  maps  the  presence  of  a 
structure  extended  to  the  north,  which  is  also  evident  in  the  ve¬ 
locity  map  (Fig.  |A.6|y),  at  slightly  red- shifted  velocities  of  50- 
80  km  s-1.  The  narrow  core  of  the  CO(l-O)  emission  in  IRAS 
23365+3604  traces  a  disk  of  molecular  gas  with  a  major  rotation 
axis  oriented  at  a  PA  ^  -55  deg. 


Appendix  A.4:  Mrk273 

Mrk  273  is  a  Sy2-ULIRG  with  a  double  nucleus  in  near-IR  im¬ 
ages  and  a  44  kpc  long  tidal  tail,  which  are  signatures  of  merging. 
The  physical  conditions  of  the  molecular  gas  in  this  galaxy  sug¬ 
gest  the  presence  of  prodigious  starburst  activity,  fed  by  large 
amounts  of  cold  molecular  gas  ([Fapadopoulos  et  ah]  [201 2a). 
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Fig.  A.7 :  Left:  IRAM  PdBI  cleaned  map  of  the  CO(l-O)  blue  and 
red-shifted  wings  of  Mrk  273,  integrated  within  ve  (-800,  -400) 
km  s-1  and  ve  (400,  900)  km  s-1  and  combined  together.  Con¬ 
tour  levels  correspond  to  lcr  (lcr  rms  level  is  0.2  mJy  beam-1). 
Right:  Real  part  of  uv  visibilities  as  a  function  of  the  uv  ra¬ 
dius  for  the  combined  CO(l-O)  wings,  with  visibilities  binned 
in  baseline  steps  of  25m.  The  red  curve  is  the  best  fit  with  a 
circular  Gaussian  model. 


by  our  new  PdBI  observations.  However,  we  emphasise  that  in 
our  analysis  we  have  chosen  to  be  conservative  and  to  define  the 
red-shifted  wing  only  from  v  >  400  km  s-1. 

As  for  IRAS  23365+3604,  we  estimate  the  size  of  the  out¬ 
flow  using  the  uv  plot  of  the  two  wings  combined  together,  which 
is  shown  in  the  right  panel  of  Fig.  |A.7|  The  partially  decreasing 
trend  of  the  uv  plot  is  inconsistent  with  a  point  source,  indicat¬ 
ing  that  the  outflow  is  marginally  resolved,  despite  the  low  spa¬ 
tial  resolution  of  our  D-conf  data.  The  fit  with  a  simple  circular 
Gaussian  (red  curve  in  the  plot)  provides  an  approximate  size  of 
(FWHM)of  1.1  ±  0.5  kpc. 

For  comparison,  we  show  in  Fig.|A.8|the  emission  from  the 
narrow  core  of  the  CO(l-O)  line,  integrated  within  ±250  km  s-1. 
This  component  traces  a  central  concentration  of  gas,  which  is 
in  regular  rotation,  with  a  major  axis  oriented  with  a  PA  of  about 
-75  deg.  (panel  b  of  Fig.  |A.8|).  The  uv  plot  of  the  narrow 
core  emission  can  be  modelled  with  a  circular  Gaussian  with 
a  FWHM  of  2.026  ±  0.046  kpc  (panel  (c)  of  Fig.  |A^). 


Fig.  A. 8:  Panel  (a):  IRAM  PdBI  continuum- subtracted  map  of 
the  CO(l-O)  narrow  core  emission  in  Mrk  273.  Contours  corre¬ 
spond  to  30cr  (lcr  rms  level  is  0.2  mJy  beam-1).  The  inset  (b) 
shows  the  first  moment  map  obtained  within  ve(-100,  100)  km 
s-1  by  applying  a  flux  density  threshold  of  8  mJy.  The  size  of  the 
inset  is  14  x  14  arcsec.  Panel  (c):  uv  plot  for  the  narrow  core, 
with  visibilities  binned  in  uv  radius  steps  of  10m,  fitted  with  a 
circular  Gaussian  model. 


The  ionised  and  neutral  gas  kinematics  is  complex  in  this  ob¬ 
ject,  characterised  by  kpc- scale  outflows.  In  particular,  Rupke  & 
Veilleux  (  2013b)  detect  two  super-bubbles  of  outflowing  ionised 


gas,  aligned  in  the  north-south  direction,  which,  due  to  the  large 
velocities,  are  likely  generated  by  an  AGN. 

The  map  of  the  blue  and  red- shifted  CO(l-O)  wings,  inte¬ 
grated  respectively  within  the  velocity  ranges  ve(-800,  -400)  km 
s-1  and  ve(400,  900)  km  s-1  and  combined  together,  is  reported 
in  the  left  panel  of  Fig.  |A.7|  We  note  that  a  significant  fraction 
of  the  outflow  (especially  the  red- shifted  wing,  see  also  Fig.  |4fi), 
is  elongated  towards  north.  Notably,  despite  the  narrow  band¬ 
width  available  at  that  time,  the  CO(l-O)  red-shifted  component, 
peaking  5  arcsec  to  the  north,  had  already  been  observed  by 
|Downes  &  Solomon  (1998),  although  only  up  to  velocities  of 
~350  km  s-1.  These  authors  stressed  that  the  velocity  and  orien¬ 
tation  of  this  structure  are  inconsistent  with  the  nuclear  disk  and 
its  east-west  velocity  gradient.  This  consideration  by  |Downes| 
&  Solomon  (|1998|)  supports  our  hypothesis  (also  discussed  in 


Section  3)  that  the  asymmetry  of  the  CO(l-O)  line  at  ve(200, 
400)  km  s-1  traces  outflowing  gas.  Indeed,  the  data  reported  by 
these  authors  only  recovered  velocities  v  <  350km  s-1,  and,  con¬ 
sequently,  they  could  not  be  observing  the  high  velocity  wings 
at  v  >  400km  s-1  which  are  instead  revealed,  for  the  first  time, 


Appendix  A.5:  IRAS  F23060+0505 


Fig.  A.9:  Left:  IRAM  PdBI  cleaned  map  of  the  CO(l-O)  blue- 
shifted  and  red-shifted  wings  of  IRAS  F23060+0505,  integrated 
in  the  velocity  ranges  ve  (-500,  -300)  km  s-1  and  ve  (300,  1 100) 
km  s-1,  and  combined  together.  Contours  correspond  to  lcr  ( 
lcr  rms  noise  is  0.2  mJy  beam-1).  Right:  Real  part  of  visibili¬ 
ties  plotted  as  a  function  of  the  uv  radius  for  the  CO(l-O)  broad 
wings.  Visibilities  are  binned  in  baseline  steps  of  30m.  The  red 
curve  is  the  result  of  the  best  fit  with  a  circular  Gaussian  model. 


Fig.  A.  10:  Panel  (a):  IRAM  PdBI  continuum- subtracted  map 
of  the  CO(l-O)  narrow  core  emission  in  IRAS  F23060+0505. 
Contours  correspond  to  5<x  (lcr  rms  level  is  0.3  mJy  beam-1). 
Panel  (b)  shows  the  first  moment  map  obtained  within  ve(-200, 
200)  km  s-1  by  applying  a  flux  density  threshold  of  8  mJy.  The 
size  of  the  inset  is  9  x  10  arcsec.  Panel  (c):  uv  plot  for  the 
narrow  core,  with  visibilities  binned  in  uv  radius  steps  of  10m, 
fitted  with  a  circular  Gaussian  model. 
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The  ULIRG  IRAS  F23060+0505  hosts  a  heavily  obscured 
AGN  (type  2  QSO).  Previous  studies  of  the  morphology  and 
kinematics  of  its  extended  narrow  line  region  indicated  the  pres¬ 
ence  of  outflows  of  ionised  gas  through  blue- shifted  asymmetries 
of  the  emission  line  profiles  (jWilman  et  al.|1999~|. 


In  IRAS  F23060+0505,  we  define  the  blue  and  red- shifted 
wings  of  the  CO(l-O)  line  within  the  velocity  ranges  (-500,  - 
300)  km  s-1  and  (300,  1100)  km  s-1;  because  the  red  wing  is 
detected  at  a  significance  below  8 cr  (Fig.  [5jc-d)),  we  combine 
them  together  to  estimate  their  size;  the  map  of  the  two  wings 
merged  together  and  the  corresponding  uv  plot  are  shown  in  Fig. 
|A.9|  From  the  latter,  by  performing  a  fit  with  a  simple  circular 
Gaussian  function,  we  estimate  their  spatial  extension  to  be  8.1 
±  2.9  kpc  (FWHM). 

The  cleaned  map,  the  first  momentum  map  and  the  corre¬ 
sponding  uv  plot  for  the  narrow  core  of  the  CO(l-O)  emission,  in- 

a-c). 


A.10 


tegrated  within  (-250,  250)  km  s  ,  are  shown  in  Fig 
The  inferred  flux  of  this  CO(l-O)  narrow  component  is  Sm  corf 

shows 


(12.06  ±  0.85)  Jy  km  s  .  The  inset  (b)  of  Fig. 


A.10 


that  the  narrow  core  of  the  CO(l-O)  emission  traces  molecular 
gas  in  regular  rotation  in  a  disk,  with  a  major  axis  of  rotation 
oriented  with  a  PA  ~  0  deg.  A  fit  to  the  uv  plot  with  a  circular 
Gaussian  function  provides  an  estimate  of  the  size  of  this  nar¬ 
row  component  of  about  6.96  ±  0.52,  which,  within  the  errors, 
is  comparable  to  the  extension  of  the  broad  CO  wings  in  this 
source. 


Appendix  A. 6:  Mrk876  (PG  1613+658) 


Fig.  A.  1 1 :  Left:  IRAM  PdBI  cleaned  map  of  the  CO(l-O)  blue 
and  red- shifted  wings  of  Mrk  876,  integrated  in  the  velocity 
ranges  ve  (-500,  -300)  km  s-1  and  ve  (400,  1700)  km  s-1,  and 
merged  together.  Contour  levels  correspond  to  lcr  (lcr  rms  level 
is  0.2  mJy  beam-1).  Right :  Real  part  of  uv  visibilities  as  a  func¬ 
tion  of  the  uv  radius  for  the  combined  CO(l-O)  wings.  Visibil¬ 
ities  are  binned  in  baseline  steps  of  30m.  The  red  curve  is  the 
result  of  the  best  fit  with  a  circular  Gaussian  model. 


Despite  being  classified  as  a  LIRG,  the  total  IR  luminosity 
(log(LIR)  =  11.97,  Veilleux  et  al.  2009|)  of  the  powerful  QSO 
Mrk  876  (PG  1613+658)  makes  it  very  close  to  being  a  ULIRG. 
The  detection  of  a  molecular  outflow  is  unclear  in  Mrk  876,  and 
the  quality  of  our  data  only  allows  us  to  set  an  upper  limit  on  the 
mass  and  outflow  rate  in  this  source.  Although  the  PV  diagram 
along  the  major  rotation  axis  of  this  source  does  not  present  any 
clear  evidence  for  high-velocity  outflowing  gas  (Fig.  the 
map  of  the  blue  and  red-shifted  CO  wings  combined  together, 
reported  in  Fig.  |A.11|  shows  some  emission  detected  at  10cr.  A 
Gaussian  fit  to  the  uv  visibilities,  binned  in  steps  of  30m,  pro¬ 
vides  a  size  (FWHM)  of  this  broad  CO  component  of  7.1  ±  1.4 
kpc. 
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Fig.  A.  12:  Panel  (a):  IRAM  PdBI  continuum-subtracted  map  of 
the  CO(l-O)  narrow  core  emission  in  Mrk  876.  Contours  corre¬ 
spond  to  5<x  (lcr  rms  level  is  0.2  mJy  beam-1).  Panel  (b)  shows 
the  first  moment  map  obtained  within  ve(-200,  200)  km  s-1  The 
size  of  the  inset  is  11  x  10  arcsec.  Panel  (c):  uv  plot  for  the 
narrow  core,  with  visibilities  binned  in  uv  radius  steps  of  15m, 
fitted  with  a  circular  Gaussian  model. 


The  narrow  core  of  the  CO(l-O)  emission,  although  detected 
at  high  significance,  is  not  completely  resolved  by  our  observa¬ 
tions,  as  shown  by  the  cleaned  map  and  corresponding  uv  plot 
shown  in  Fig.  |A.12|  The  inset  (b)  of  Fig.  |A.12  shows  that  the 
central  concentration  of  molecular  gas  in  Mrk  876,  responsible 
for  the  bulk  of  its  emission,  is  likely  rotating  in  a  regular  disk 
or  ring,  with  a  rotation  major  axis  oriented  at  a  PA  =  -37  deg. 
We  infer  for  this  component,  within  (-250,  300)  km  s-1  (i.e.  the 
maximum  velocity  extent  of  the  rotation  as  evidenced  by  the  PV 
diagram  in  Fig.  [6J:),  a  physical  size  of  7.18  ±  0.46  kpc,  which  is 
perfectly  consistent  with  our  estimate  of  the  size  of  the  putative 
broad  CO  wings. 


Appendix  A.7:  iZw  1  (PG  0050+124) 


Fig.  A.  13:  Panel  (a):  IRAM  PdBI  continuum-subtracted  map  of 
the  CO(l-O)  core  emission  in  I  Zw  1.  Contours  correspond  to 
10cr  (lcr  rms  level  is  0.2  mJy  beam-1).  Panel  (b)  shows  the  first 
moment  map  obtained  within  ve(-200,  200)  km  s-1,  by  applying 
a  flux  threshold  of  3  mJy.  The  size  of  the  inset  is  16  x  16  arcsec. 
Panel  (c):  uv  plot  for  the  narrow  core,  with  visibilities  binned  in 
uv  radius  steps  of  5m,  fitted  with  a  power  law  function  model  of 
the  form  oc  r-2. 


In  the  LIRG  and  luminous  QSO  host  I  Zw  1  (PG  0050+124) 
we  do  not  detect  any  CO(l-O)  broad  emission  component,  de¬ 
spite  the  fact  that  its  UV  and  X-ray  spectra  show  some  evidence 
for  high  velocity  (~2000  km  s-1)  ionised  outflows  (Costantini 


|et  al.|2007).  Figure  [A 713]  shows  the  map  and  uv  plot  of  the  nar¬ 
row  CO(l-O)  line  of  I  Zw  1,  integrated  within  v  =  +  250  km  s-1. 
The  CO  line  emission  appears  well  resolved  by  our  observations, 
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and  the  uv  visibilities,  binned  in  baseline  steps  of  5m,  can  be  well 
fitted  with  a  power  law  model  of  the  form  oc  r-2,  which  gives  a 
FWHM  of  1.258  ±  0.041  kpc.  The  molecular  gas  traced  by  this 
narrow  CO(l-O)  component  is  mostly  rotating  in  a  central  disk 
or  ring,  with  a  major  axis  of  rotation  oriented  with  a  PA  =  -50 
deg  (inset  (b)  of  Fig.  |A.13|). 
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